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ABSTRACT 

The Sand river map Sheet, located in east-central 
Alberta, possesses several bedrock channels buried beneath 
Quaternary glacial deposits and known only from drill hole 
information. A total of 995 samples from 78 drill holes were 
analyzed for 22 compositional properties in an attempt to 
delineate the stratigraphy of the till units present. The 
measurements inciuded granulometry, grain lithology, calciun 
carbonate equivalent and bulk (major and minor element) 
chemistry. A reliable, rapid, constant -volume carbonate 
analysis apparatus was developed for the project. Bulk 
Chemistry was performed by fusing till samples into glasses 
and analyzing them using energy-dispersive electron 
microprobe techniques. The correction procedures required 
for quantitative energy-dispersive analysis are described in 
some detail, along with EDATA, a computer program for 
applying then. Corrections are applied for atomic number, 
absorption and fluorescence effects, instrumental 
calibration corrections, background (X-ray continuum) 
Shaping and fitting as well as peak overlaps. 

Problems are often encountered in contour mapping of 
irregularly spaced data points by computer, and new 
techniques has been devised and implemented in the TRIMAP 
package, based on subdividing the map area into triangles. 
Features include the use of triangular element data 


Structures, a local homogeneous coordinate system, the 
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automatic generation of a triangular network, the selection 
of a suitable triangle optimization criterion and the choice 
of a suitable smooth interpolant over individual triangular 
domains. 

Four distinct till units were distinguished in the 
study area on the basis of their compositions. A high 
carbonate middle till separates an upper and a two member 
lower tili. The compositional parameters of these tills were 
categorized as textural, erratic or local. In the fine 
fraction of a tiil unit, vertical compositional variation 
was low by comparison with lateral variation, which appeared 
to be related to sub-ice topography. The glacial history 
proposed relates the four suggested ice advances to the 
varying regional drainage as indicated by the major bedrock 
channels. The regional topography at each stage was 
estimated, Compositional analysis of tills, as weli as the 
topographic reconstruction of stratigraphic boundaries, is 


important in studies of glacial history. 
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41.0 INTRODUCTION 

For many years the glacial geology of the Prairies of 
Western Canada has been known primarily in two dimensions. 
This was because of the relatively infrequent occurrence of 
Suitable vertical sections. Consequently, although the 
glacial stratigraphy over much of the Prairies might be 
expected to be relatively consistent and extensive, little 
regional information has been available below the uppermost 
units. 

A general exception to this state of affairs has been 
the work of Christiansen (1968a,b, 1972a,b) utilizing 
information and samples from holes drilled in the process of 
regional hydrogeological surveys. He was able to define the 
regional glacial stratigraphy over much of central 
Saskatchewan by defining the stratigraphic relationship 
between drill holes on the basis of electric logs and sample 
examination. Sauer (1974) mentioned the practical value of 
till stratigraphy in highway engineering and water supply 
problems. In particular, where the tills differed in their 
engineering properties, embankment failure could have been 
caused by water seepage along inter-till boundaries, and 
major surficial aquifers were usually located between or 
below major till layers. 

Visual examination of available vertical sections is of 
great value in conventional field work, but laboratory 
analyses are often very helpful in the identification of 


till units. When field examination is restricted to samples 
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taken from a hole during drilling operations, laboratory 
analysis is necessary to determine criteria for 
distinguishing superimposed till units. 

The primary objective of this study was to examine the 
regional glaciai stratigraphy in an area of Alberta where 
drill hole samples and data were readily available. The area 
was also of interest as three bedrock channels were known to 
be present, and it was hoped that some relationship could be 
observed between regional glaciation and drainage. It was 
therefore necessary to determine those compositional 
properties of greatest value in till identification as well 
as to determine efficient sampling and analytic strategies 
for any future work. In the pursuit of these objectives it 
became necessary to develop new analytic tools. The 
development of improved apparatus for measuring till 
carbonate content was one result (see Chapter 4), as was the 
development, in association with Dr. D. G. W. Smith, of 
fully quantitative energy dispersive microprobe correction 
procedures as used in electron microprobe analysis of the 
bulk chemistry of the fine fraction of till samples (Chapter 
5). Finally, in the display of topographic data, grid-type 
computer contouring algorithms were found to be inadequate 
because of their inability to resolve features with few 
sample points, and the difficulty in obtaining a surface 
that passes precisely through all the data points. A new 
approach to contour mapping was consequently developed that 


honoured all data points instead of smoothing the surface, 
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as well as permitting the user-definition of features such 
as valley-lines. 

Due to the variety of the subject matter covered, some 
explanation of the order of presentation may be applicable. 
After a brief description of the study area, the laboratory 
methods used in the examination of the drill hole samples 
are discussed, along with the statistical tests used to 
evaluate the various analyses. Because the carbonate and 
microprobe analysis techniques were developed especially for 
this study,the next two chapters discuss them in detail. In 
a Similar manner, the display and interpretation of the 
greater than 20,000 analytic results required the 
availability of suitable display methods, and hence the 
chapter on contour mapping precedes the detailed examination 
of laboratory results. Finaliy, one of the primary purposes 
of the laboratory analyses is the delineation of the till 
units found. Thus an interpretation of the glacial history 
of the area reguires in advance both a knowledge of the till 
units present, as well as some background knowledge of the 
preglacial history. Historical geology thus represents the 


Synthesis of much previous work. 
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In 1974 the Groundwater Branch of Alberta Environment 
published a report entitled "Buried Channels in the Edmonton 
- Lac La Biche - Cold Lake Area, Alberta," (Yoon and Vander 
Pluym, 1974). This was a geological study of three well 
defined buried bedrock channels: the Beverly, the Kikino and 
the Helina. Field work in the area included test drilling, 
lithologic and eiectric logging, and side-wall sampling of 
the drilled wells. One hundred and fifty two test holes were 
drilled, usually to bedrock, with sidewall samples taken 
from ninety eight of them. Additional well log information 
was obtained from several other sources, including oil and 
gas exploration wells. 

In this study data were utilized from seventy eight of 
these drill holes in an area of good sampie density and 
potentially interesting geology where three buried channels 
were present. All samples from these drill holes were 
examined and 995 were analyzed in the iaboratory. Fig. 1 
shows the location of the study area and its reiation to the 
previous report. All drilled holes were located between 
townships 56 and 68, ranges 2 to 17 west of the fourth 
meridian. The study BG Relic mostly within the Sand River 
(73L) 12:250,000 map sheet of the National Topographic 
Survey, although a few drill holes are located just to the 


south {map 73E) or west (831) of this area. 
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Fig. 1 Index map of the study area. 
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2:2 Topography and Drainage 

The study area is part of the plains region of east- 
central Alberta. Topography generaliy rises to 700 metres in 
the southwest and slopes gently to the north and east. 
Moderately higher ground is encountered in the northeast 
corner of the map area. The low points (below 550 m) in the 
area are around Lac La Biche and between Moose and Cold 
Lakes. A broad moderately low region connects them (Fig. 2). 
Three major drainage systems are present in the area, with 
little in the way of high ground to separate them. The La 
Biche River drains the low ground in the northwest and is 
part of the Athabasca and Mackenzie system, flowing into the 
Beaufort Sea. The Sand and Beaver Rivers drain to the east 
via Coid Lake and the Churchill River into Hudson Bay at 
Churchill. The North Saskatchewan River follows a more 
southerly route and drains into Hudson Bay at York Factory. 
The area is covered predominantly with ground moraine 
(having a local relief of less than 3 m) or hummocky moraine 


possessing more pronounced mounds and depressions. 


revious work in the area includes the report of Yoon 
and Vander Pluym (1974) on the buried channels, along with 
studies of the bedrock topography by Carlson (1967), Lennox 
and Carison (1967) and Carlson and Currie (1974). Borneuf 
(1973) compiled a hydrogeological map of the Tawatinaw map 


area, which includes the western part of the area of this 
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study. Le Breton (1963) did the same for east-central 
Alberta. Bedrock geology has been described by Carrigy 
(1972) and Green (1972). A hydrogeological map of the map 
Sheet immediately to the east has been produced by 


Christiansen and Whitaker (1974). 
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3-1 Types of Analyses 

Four categories of anaiysis were performed on each of 
995 samples from the map area: grain size analysis, 
lithological analysis by grain counting, carbonate analysis 
and bulk (major and minor element) chemistry. Complete 
details of the laboratory procedures used are found in 
Appendix 1. A general discussion of each is given below 
along with a statistical summary. Table 1 gives the legal 
locations of the 78 drili holes and Appendices 2 and 3 give 
the results of the analyses performed for all 995 samples. 
Appendix 4 plots these results beside the electric log for 


each drill hole for ease of examination. 


The samples used in this project were obtained from 
holes drilled, logged and sampled by Alberta Environment 
staff. While it would have been desirable to obtain more 
frequent and regularly spaced samples from each hole, 
drilling schedules, the inability to sidewall-sample coarse- 
grained material and the occasional cave-in prevented this. 
Because the sampler may not have penetrated the side of the 
hole precisely where it was intended it was important that 
each sample was described in the laboratory prior to 
processing. Comparison with the driller's and electric logs 


identified any mis-located samples. Because, for example, 
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some tills, lake deposits and shales may be composed of very 
Similar materials it was of considerable importance that 
those structural features that could distinguish them were 
described prior to sample crushing. Only those samples 
considered to be till were analyzed further. 

The term till may have at least two meanings. In the 
broad sense it includes any material deposited by glaciai 
ice, and in this context extensive thicknesses of bedrock 
material that have been interpreted as having been moved by 
the ice have been classified in this study as part of a 
particular tili unit. Dreimanis (1976) discusses deformed 
bedrock and block inclusions in tili. In a more limited 
sense, till is considered to be a material posessing a 
Ccharacteristically wide range of grain sizes that has been 
deposited directly from glacial ice. Goldthwait (1971) 
discusses till as a fairly homogeneous material of mixed 
particle sizes, It was in this sense that individual driil 
hole samples were classified as till. While the first of 
these two properties was readily observable ina small 
sample, the requirement that it be ice deposited was less 
readily ascertained than would have been the case for 
exposed geological sections. Consequently there is a very 
real probability that some other diamicton was mis- 
classified as till, the most likely candidate being till 
material that was re-mobilized in some manner, and hence was 
no longer strictly till. However, this problem has been 


reduced to some extent by the scale chosen for the current 
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study. Unlike the examination of a single exposed section, 
lithologic units much less than 5m thick are often difficult 
to resolve, due both to the vertical sample spacing and the 
distance between drill holes. Thus thin diamictons 
interspersed with water-laid deposits would usually be 
assigned to a "mixed" category in the log of a drill hole, a 
not unreasonable classification for a regional-scale study. 
Thus only extensive thick diamictons would be likely to be 
stratigraphically considered tills, an interpretation that 


would probably be correct. 
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Laboratory methods were evaluated by analyzing a set of 
arbitrarily chosen test samples. The results of these, 
analyses are given in Table 2. For the grain size, grain 
count and carbonate analyses, measurements were taken 
several times on each of several samples, using different 
sub-samples of each test sample to produce the replicate 
analyses of Table 2. In the case of bulk chemistry, an 
additional step - the preparation of a fused glass of the 
Sample - was needed. Several glasses were made from each 
test sample and bulk chemical analysis was performed twice 
on each glass using the electron microprobe. May and 
Dreimanis (1976) suggest a variability model in which an 
observation is the sum of the population mean and the 
natural, sampling, preparation, analytical and random 


variabilities. Variation between the different analyses for 


q.2 


aces Sea iidodyte eet é 
{uote he) eda gis Aoras #2 


ihe Bad meee shone ee bE oF faa 


aodr her Acad anagy eztaag O20 


od) yliades Sino 2 Ake gse akan > 


a gee: LDEES #30 Re ait a2 ROR sree bis / ae 


ausiie sheer é sokpat F 40% ta eS a Hae ore 


Pa A inal vi sete ae pias ae ie ~ an seyias i B 4 rat ii 
fi 
So aebreho: bi 


( ue mes 


grad HOeTETSATISIAE 16) 42 3 Gist 
,sSaa ob abet 


' i a 3 


ahodred pyeewided 29 Poets 
Pie | bidet Yt a 
my ey aul 

7 


At is 
its es if 
“-— 14), ‘ 
| 7 


= _ 
~ } 
S2a2b yvesis Tew) pa- TO). ay 9 LS ibe osben ie 


4 


So tee so pHELELBe Yo bet aliteys 
- | fe 
opiedg 30, 27lse? SAP wae 4 abe Supt 
vom sy 
Jotay stay es aces ssn 285 ahh ‘ Fn 


ae dues Rebioaien 32 Laie y eaves a ie aA 
Ou geeotdiqeavet? soubdxg of Finate i Jeet dobar 
oa! \eeFeeaeno. Ahad ta SEM aie, te i 

ed¢ Yo .2e61y beant = 20 sonerisEsis ‘aa =< 

fpes Worst fr Bit: eta asanele taitey oe 


OLL¥ 3 ad a a atey sans {xabgeeo ae 


rma) 


De oe baa tox jeseagototm domes: sid Buleo anes yen 
7 iv > Pc) Uns mS 
ic. Gan i dazay, at fe hor (i tfedsin a veogeud (overa 2 giro : 


. : a lig hy Dis ne ites abitstoneg a3. ae aus bits, ae ¥ bit — 
ae i: | 


\ a wha Dee. pinat : nr rs 


V2 


3 “soueTeatTnbe eyeuogied antoTeo. pue AboTouygtyt uteszb 
eZTS uTeib ‘“Xijstwayo y {nq 30; seinpssoid Aioj,ero0gey 
4$93 ©} pesn setTdues 4saQ3 Jo seskTeue Jo s3ytnsay Zz STqeL 


*syabg 
L9°P 0°s 0°S $°9 9°Y sy gy ov £°Y $°2 9:2 2. ¥22 £°2 fo9e9 
a a a a ee eae ee ee 
A v ¢ z t v ¢ z T y ¢ z TT  83edT day 
¥ bs 1 1 1 x x x N ry c fg c atdsrs 
$6°<T S°LT — -0°9tee St S*1t= Orc St OAT — Ont OLE PAE "210 
$0°C $°z $2 s‘y Oy . OF $°0 o'T oT $70 (O'T qty 
00°s st $°0 s°t o't ot (hay 0°2 oz 05¢2.. SES dIStA 
0¢°69 o°0L O°SZENS°69 mS: 20ers 07 sso §°S9— (09 S°7Z 9°69 Ty 
SRST $8 OT iA ae OEE Ges 426%.) S30¢ 0°6Y Ost =: OT TE397 
9 
A s ? ¢ z T $ 9 £ z T s2eottdoy 
v I I I I I i Hu HN u H oT dues 
trot ; TL°2T) «(92ST £6°G)~=—GT'OL SOSZ'8 02°24 «= G8*k sé OS "ck OLecT Anta 
95° tS O9°FS -GS'SSe7-StocSin efcSS. lcclSa feb SESS Get. — 7550S — 16s VIS 
L8°S 35,9) mpelc9 me ch-O em St: yee 970 19°07 se*§ ED 98°» 69°S  Ofc-O¥T 
8c°ST Spot LL°IT ECT «OOTIT «Ob PE «= TNEZ——sHOTIT. «COPE. = «OT «E9°ST = OF T-09 
60°6 Zorba OY Limerce Olan tecO Camel SO G0 -Smm Ole = S6a6, €T°Ol = $8°Z 09-09 
wi°S z8*y = %S°e = TTS) =—o9S*t C'S ess ao"k 0 69°S ORF OF-92 
ster S2°f) WEST es BCS ae SiS eee SS -¢: L9°b -SS*b | Oke $0°9 svt" 0z-O1 
9 : 
A z T z T z T z T z TL a3eattday 
Y 9 u) d d | a a a a) 2) atercs 
a a a ee I yy ee Oe ey eee 
-SO°F SETH SOF BSZ Pr B90'R 98S LZ0°t O8T°h B8I'y 90T'h SLRS ScG’d If0°h LO6°S eFL°S CABS GESTS aj 
tore USED FESO | EERO CGEE'O: ©«OOGP"D = S7S°O—s BGO = 00S°0 Ss 6Z¥°O SS OOBr°OEO )8=8s«#9SS°0 = 9r°0 «= dzP0 =O 0S°0)— SOF tL 
S505. <tOyl: ACOOTS = ~<GT°S GSE7T. 906° er0-% «.69D%S G10“, SNS. “GDC: 2616 T «= «1E27Z 4 G56 OSS CSc | SZ1sC $9) 
CuST ELOY WOU cel «OBL «GEST «OSS8*E «= SOG*E = WEO°Z «EEG. SZG°Y = BOGE «= ssAZG"I. «OC SSB*T «= GG"T 989°) 6BB*T a 
TSB°LE LEGO"SE LIS EL FOES LCS LE Clete TSL°SE CLLEE CLS*SE TER°EL ISO'HE STZ we CLOSE BLT VE SLE°rS WLS teP'EE TS 
SOLE GORSoe: ZS8°E |  BES*L BOLL OSL ALE Eek VEER SL GUE 85902 L090 ee ATS StS*2  SZOTL  TSEtLe TOLRL Vv 
sec" £T2°G  «btL°OGEL°OTOB'O969°0)SFL4L"0 GLZ°0 ss £69°0 « LTL"O so LPSO) SOOELZ*O)=—«(9GZ°0 Ss 299°0 Ss OF9°0SsSTB*O OC 9"0 Bx 
Wl 055°0 070 3=—s- $970 0°0 0°0 TOr°0 SzT°O «Reza ~—«BS0°0 0°0 00) = SET" 0°0 0°0 oo 992°0 tN 
9 < T 2 t z t Fay aes 2 7 z T z T z T Sys{{ouy 
4 ¢ ¥ ¢ £ z z t t ? ¥ £ £ z rd t t Sst79 


v y q g q & @ a @ y Vv y Vv vy v ¥v v a{ars 


6a. 72 
zlesds> 


eanzp +773 


“-_ goaeierivps 9¢sa© 


——————— 
62 enacey 
Mae) 2278 vos. 


sues 


z ae 
ons SMO, ALC 
ths wetl, SHE . 

Rags 


Wa etor 


¢.> 


dxys> sviciso 


eae: “abt, 
D Ee Re + Eo ae ; = 


ng, le’ @ 


cag) fe Tet 


ae ee 


eat te ee 


the bulk chemical analysis may therefore be due to real 
variation between the samples analyzed (natural plus 
sampling variability, or between-sample variance in this 
example), to compositional variation between material taken 
from the sample to make each glass (preparation variability 
Or sub-sample variance) or to variation between individual 
measurements on each glass (analytical plus randon 
variation, or error variance). It should be noted that it 
was not known in advance whether a particular set of samples 
could be expected to differ significantly in all the 
analyses performed on it, although a different choice of 
samples might well have shown such a difference. Thus the 
ability to distinguish samples in the test set is of less 
interest than the standard deviation of the within-samples 
(or between glasses, within-samples) variance. Analysis of 
variance techniques were used to examine the various 
components of the variation in laboratory analyses on the 
test samples.In the case of the bulk chemistry, two 
arbitrary samples, A and B, each had 4 glasses (sub-samples) 
made from them. Two analyses were performed on each glass. 
The analysis of variance model is a one way classification 
with equal sub-sample numbers (Steel and Torrie, 1960, 
Section 7.11). The analyses are given in Table 3. The within 
sample variance may be ascribed to both variation between 
glasses made from the same sample and variation in multiple 
analyses of each glass; hence two F tests may be made, 


comparing variation between samples with variation between 
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glasses, and comparing between glasses and within glasses 
variance. In the second case significant differences between 
glasses for the same sample were not found, showing that 
glass preparation was not a significant source of error. 

Table 4 gives the analysis of variance for the grain 
size, grain count and carbonate analysis. The model is a 
one-way classification with equal replication (Steel and 
performed on each sub-sample, the within-sample variance is 
a combination of sampling and error variance. Duplicate 
analyses are not performed on each sub-sample in carbonate 
analysis because the method destroys the material used. For 
grain counts, analytic error is only found if the operator 
classifies a grain differently on two different passes, a 
relatively infrequent occurrence. In addition, the small 
humber of suitable 1-2mm sand grains present in a sample 
necessitate recombining the analyzed sub-sample with the 
test sample prior to selecting a further 200 grains for 
analysis. While it may be possible to perform grain size 
analysis several times on each sub-sample this was judged to 
be technically awkward, especially for the silt and clay 
fractions. 

Details of which particular compositional variables 
differ significantly between test samples are given in the 


relevant sections below. 
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Material from the lightly crushed sample that did not 
pass through a U.S. standard #230 sieve (0.063 mm) by dry- 
sieving was washed, dried and sieved through #10 (2 mm), #20 
(1 mm), #40 (0.5 mm), #60 (0.25 mm), #140 (0.125 mm) and 
#230 size sieves. Dry material initially passing through the 
#230 sieve therefore consisted of silt and clay, and pipette 
analysis was used to determine the clay content (see Day, 
1965). The pipette method was chosen over the hydrometer 
method due to the reduced labour invoived. In the former 
once particle settling has commenced a single measurement is 
reguired after the calculated settling time, whereas in the 
latter method readings must be taken at systematic 
intervals, severely limiting the number of samples that may 
be analyzed at one time. In addition, the pipette method is 
more readily usable on smaller samples than the hydrometer 
method. Material between 2mm and 0.063 mm was considered to 
be sand, between 0.063 and 0.002mm silt and below 0.002mm 
clay. 

In the analysis of variance tests described earlier the 
test samples could be distinguished by the weight 
percentages of the 1-2 mm. size fraction at the 95% level, 
the 0.5 - 1 mm. fraction at the 90% level and the clay 
fraction at the 95% level. The intermediate size ranges did 
not differ significantly between the test samples. Although 
well able to distinguish the test samples, the absolute clay 


percentages must be suspect due to the long storage times 
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and consequent disaggregation problems associated with the 
Sidewali samples used (see Appendix 1). Pawluk and Bayrock 
(1969), for example, found an average of 30% clay (less than 
2 Micron material) for 20 samples from widely spaced 
locations in the Sand River map sheet, markedly higher than 
the 10% observed in the present work. Nevertheless the 
coarse sand and clay fractions appear to be of use in sample 
differentiation due to their significant F tests mentioned 
above. It should be noted that the non-standard 
granulometric methods descriked above were implemented due 
to both the large number and small sizes of the available 


samples. 


Lithological analysis is commonly used in attempts to 
distinguish between till sheets.Some results from the 
Edmonton area are quoted in Westgate et al. (1976). His 
results were obtained by counting boulders in the field or 
pebbles in the laboratory. Because the size of sidewall 
samples precluded the use of these size fractions, the 1 - 2 
mm sand fraction was used instead. Dreimanis and Vagners 
(1969,1971) discuss the use of differing size fractions for 
examining the variation in lithology in a till, and Willman 
et al. (1963) use the mineralogy of the sand sizes to 
distinguish between till sheets. They concluded that sand 
lithology may successfully distinguish till sheets, but care 


must be taken not to confuse measurements obtained in one 
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area using one size fraction with results obtained elsewhere 
using a different size fraction. 

The method used involved the classification of 200 
grains from each sample into five categories: local shales 
and sandstones, acid igneous rocks, basic igneous rocks, 
carbonate rocks and undifferentiated quartz grains. Five to 
Six samples could be analyzed per hour if interspersed with 
work less tiring to the eyes. The test samples could be 
distinguished at the 95% confidence level by the local and 
carbonate rocks and at the 90% level by the acid and basic 
igneous rocks (Table 4). AS shown by the fairly high 
coefficient of variation for the basic and carbonate grains, 
the number of counts observed were rather low for 
assumptions about the normal distribution of the results to 
be certain. Statistics based on these two measurements 


should therefore be treated with caution. 


Among the rock types that were probably traversed by 
the glacial ice are the limestone and dolomite found in the 
belt of Devonian rocks fringing the Precambrian shield to 
the north and east of the map area (Douglas, 1968 and Green, 
1972). Rather than counting sand grains to determine the 
lithology of the parent rocks, however, the content of 
carbonate rocks is readily determined by mixing the till 
with hydrochloric acid and determining the amount of carbon 


dioxide evolved. Dreimanis (1962) and Dreimanis and Vagners 
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(1969, 1971) used this method to determine calcite and 
dolomite content of tills in Ontario. Christidnsen (1968a,b, 
1972a,b) used the same approach to measure the total 
carbonate content of tills in Saskatchewan. Along with other 
workers they measured the volume of gas at constant 
pressure. Dreimanis and Vagners used a manually balanced 
glass manometer, while Christiansen and Ross (1971) used 
electro-mechanical equipment of their own design. The advent 
of relatively inexpensive solid-state pressure sensors 
permitted the development of a constant-volume device that 
is inexpensive, reliable and permits rapid analysis even 
when used by unskilled operators (see Chapter 4). In some 
parts of North America limestone and dolomite can be 
distinguished on the basis of their rates of solution, 
(Dreimanis, 1962). This has also been attempted in Alberta 
(Westgate, 1972). The current work, however, suggests that 
this distinction is not easily made and like Christiansen 
(1972b), only total carbonate or total carbon dioxide 
evolved should be used, at least in central Alberta. 
Statistical tests using total carbonate gave excellent 
Sample differentiation in the test analyses. The between- 
samples variance was significantly greater than the within- 
Samples variance at better than the 99% level (Table 4). 


Further details of the equipment are given in Chapter 4. 
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Bulk chemistry of the silt and clay fraction of till is 
potentially as useful as carbonate analysis for 
stratigraphic correlation. May and Dreimanis (1973) have 
separated tills in southern Ontario on the basis of trace 
element analysis, and Pawluk and Bayrock (1969) analyzed 
Alberta tills for some trace elements, but little work has 
been done on the major and minor rock-forming elements. 
Energy dispersive analysis on an electron microprobe was 
used to analyze the till samples for 22 elements 
Simultaneously, 8 of which occurred within the detection 
jimits of the method. Details of the methods developed are 
given in Chapter 5. Sample fractions passing through a #230 
Sieve were fused into glasses using the techniques of 
Schimann and Smith (1976) and analyzed in the microprobe. 
Raw data were processed using the procedure described in 
Gold and Smith (1976). With the exception of Na, errors due 
to preparation of glasses were unimportant on the 
Statistical test data. The samples could be distinguished by 
Na at the 90% and Fe at the 95% confidence levels, but the 
Na results were considered suspect because of the high 
coefficient of variability (Table 3) and the raw analysis 
values being close to the detection limit for that element. 
Due to the closeness of their analytic peaks and the low 
intensities observed it is possibile that the Ti figures also 


include minor amounts of Ba. 
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4.0 CARBONATE APPARATUS 


The normal method of determining the carbonate content 
of a soil or till sample requires that the carbon dioxide, 
generated by reaction with dilute (20%) HCl, be measured by 
the change in volume at constant (atmospheric) pressure. 
This is usually done by manually baiancing the pressure on 
either side of a U-tube (Dreimanis, 1962). Other work, 
attempting to relate the rate of solution to the carbonate 
composition and grain size, has been done by Skinner and 
Halstead (1958), Skinner et al. (1959), Turner (1960) and 
Turner and Skinner (1960). Small negative pressures should 
always be maintained except when readings are being taken, 
in order to reduce the chances of leakage. Temperature 
corrections applied using the gas law, PV = nRT are usually 
Only applied to the generated gas volume, and ignore the 
expansion of the dead air in the system, which can cause 
errors of the same order of magnitude as the measured 
reading. This is important since the initial heat of 
reaction is dissipated through the glass walls of the flask 
in the course of the analysis. Manual reading of volumes is 
difficult, especially when they must be taken at short time 
intervals. More elaborate constant-pressure electro- 
mechanical devices have been developed (Dyck and Perkins, 


1974). 
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4.2 Equipment 

The equipment developed here detects the carbon dioxide 
evolved by monitoring the pressure change inside a system of 
constant volume. If pressures are kept low (a few p.s.i.) 
the effect of pressure variation on the reaction or the 
solubility of carbon dioxide in the acid will be minimized. 
The pressure is monitored by a small recently-developed 
solid-state detector from National Semiconductor and the 
output fed directly to any chart recorder (Fig. 3a). The 
detector's sensitivity is much greater than is required for 
this apparatus. Use of a chart recorder permits use of a 
stronger acid, usually 50%, and hence a more rapid reaction, 
easy observation of the end of the reaction and of any 
errors such as a leak in the flask seal or a slow stirrer. 
Analysis of typical till samples showed no additional gas 
evolution for acid concentrations between 20% and 75%, 
suggesting that for samples used in this study solution of 
other than carbonate compounds was not a problem. It should 
be noted that in other areas the use of high acid 
concentrations may result in the solution of till 
constituents (€.g. sulphides) not noted in this study. 

The detector is connected to an acrylic flask with a 
screw-down sealing cap and an internal acid reservoir (Fig. 
3b). A 1.00 + 0.005 gm sample and a magnetic stirrer are 
placed in the flask, the cap screwed down, the flask tilted 
to mix the acid and sample and then piaced on the stirrer. 


The screw cap minimizes the chances of leakage (which is 
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a) ELECTRICAL CIRCUITFOR CARBONATE APPARATUS 


b) SAMPLE FLASK FOR f 
CARBONATE APPARATUS | 


2 1/4" OD. 
ACRYLIC TUBING 
WITH 3/16" 
WALLS 


— 3/8" ACRYLIC 
SHEET 


Pig. 


3 Construction of carbonate analysis equipment. 


LS 


es 
or ey eis a4 
Rai SS : 
\ t 7 i 


Ae A DA Pn ry 
miem 


sian Me CL REN 
Mi daiod dd hove ft’ mys Pa ah ; STE 


Dee 5 
D a 5F in iu 
rer a 

' 0 


bt 


| Arora. 
rt gong 2 op LE CA wy 


26 


rapidly noticed if the cap is not screwed tight enough). 
With two sets of equipment and spare flasks at least eight 
samples may be run in an hour, frequently more if the 
dolomite content is low. Since the equipment is battery 
operated and portable it may be used in the field if manual 
agitation and a volt-meter are used instead of the magnetic 
stirrer and chart recorder, which both require a 110-volt 


power. 


Pyrex glass has a thermal conductivity of 26 calories 
per second per square cm for a one degree celsius per cm 
thermal gradient. In the same ane acrylic has a thermal 
conductivity of 4 - 6. Thus as the acrylic flask is a much 
better insulator than glass the heat of the reaction is 
retained in the flask with very little loss during the 
duration of the run. This is supported by the fact that, 
using the acrylic flask, no pressure drop has been noted for 
at least 30 minutes after the reaction has ceased. 
Conversely, with the acrylic flask external temperature 
changes during a run are of little concern. Magnetic stirrer 
or hand warmth is capable of penetrating glass, as was 
readily observed while monitoring pressure changes in a 
glass flask. For small changes in temperature and pressure 
the heat of reaction is approximately proportional to the 
Carbonate content of the sample, ie. 


[4.1] Te = 2, + k Moe 
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where Mp is the mass of carbon dioxide generated and Ts and 
Tj are the final and initial temperatures of the system. 


Since the volume V of the system is constant, as is M the 


a’ 
mass of dead air in the system, the ideal gas law gives the 
initial pressure of the dead air in the system as 

PPecNOnePs =OREDAtA 8A Me 
the final partial pressure due to generated carbon dioxide 
as 

P53d. HPE seROTe Mg 71.9 
and the final partial pressure due to the dead air in the 
system as 

Br dAtieh d=ek TpAn spZeue 
R being the gas constant. Substituting for Te and 
Simplifying, the final pressure difference 

(4.5] PA SOP Atepgtoee, 

= thKGteAC Dy tuk Vonguetial) joZom: 

Thus if the mass of carbon dioxide is a small fraction of 
the dead air mass the quadratic term can be ignored and 
pressure change considered proportional to carbon dioxide, 
and hence carbonate, mass. In most laboratories the 
dependence on initial atmospheric temperature and pressure 
can be adequately compensated for by regular calibration 
with reagent grade calcium carbonate. As shown in 


Christiansen and Ross (1971) and Perkins and Dyck (1973) it 


is important that stirring speed be standardized. 
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Calibration of the apparatus is direct and analyses for 


calcium carbonate equivalent may be calibrated easily by 
running a standard of known calcium carbonate composition 
every few hours, the pressure being nearly proportional to 
the carbonate content, as shown above. The standard error of 
analysis is 0.2% (see Table 4), and is partially dependent 
on the recording device used and the accuracy of sample 
weighing. The detection limit is visually estimated as 0.2%. 
Distinction of calcite from dolomite content due to the 
break in slope of their curve due to their differing 
solubilities does not appear to be feasible in the test area 
since the local limestone may take two or three minutes to 
dissolve (see Fig. 4), as opposed to the less than one 
minute quoted for limestone in other regions (Dreimanis, 
1962). While this may be due to partial dolomitization of 
the limestone samples, no break in slope was readily 
apparent in till sample analysis curves. For simplicity of 
calibration all figures are quoted as calcium carbonate 


equivalent. 
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f 3 ote TYPICAL CARBONATE ANALYSES 


. 1gm. of a typical till with tc 
- 86.3% CaCOs equivalent 


, eo een es eee ee ee ee oe 


Atmospheric 
Pressure 


100 mg. dolomitic limestone. 
{U.S. Notional Bureau of Standards 
standard sample 88a.) 


~ bbs 


Sealing 


Vn ssure 


100 mg. of local limestone 


100 mg. CaCOs reagent 


L : ee ler ass 
Fig. 4 Plots of pressure (x-axis) versus time (y-axis) 
for several typical materials, using the apparatus of 
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2-1 Principles of Electron Microprobe Analysis 
In an electron microprobe, a beam of electrons is 
focussed on a polished specimen in order to excite the atoms 

of the sample and produce the characteristic X-ray spectra 
of elements present. Ignoring various factors such as loss 
by absorption, the intensity of any particular X-ray line of 
any element is proportional to the concentration of the 


element in the sample. Thus if a standard and a specimen are 


analyzed under the same conditions: 


see i sp = Cot * Tsp/1st 
where Cp and Cot are the concentrations in the specimen and 
standard, respectively, and Isp and Ig; are the measured X- 
ray line intensities in the specimen and standard, 
respectively. 

The conditions under which a sample may be analyzed, 
however, can vary appreciably. Different instruments may use 
different angles of incidence between the electron beam and 
the sample surface, as well as different take-off angles at 
which the generated X-rays are sampled. Different operating 
conditions with the same instrument would include the 
operating voltage which accelerates the electron beam and 
the current flowing through the bean. 

Since it is necessary to measure the intensity of a 
particular X-ray line at a particular energy a technique 


Must be used to distinguish between the various lines 
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emitted. The most common approach is to diffract the X-rays 
using a suitable crystal and then rotate the detector until 
the desired wavelength is being sampled. This is described 
in Rucklidge (1976). 

A more recent approach, energy dispersive analysis, 
uses a solid-state, lithium-drifted silicon (Si(Li)) 
detector. The size of the pulses output by the detector is 
dependent on the energy of the impinging X-rays, while the 
frequency of the pulses depends on X-ray intensity (or rate 
of photon reception). The pulses are therefore catalogued 
according to their energy and the counts stored. Thus a 
typical analysis will result in a spectrum of counts in each 
of 1024 channels, each 10eV wide. While extremely convenient 
for simultaneously collecting intensities from several X-ray 
lines, the ability to resolve different energies is not as 
good as with diffraction methods. Thus characteristic peaks 
are relatively broad, causing two particular problems. One 
is that peaks may overlap each other extensively, making it 
difficult to estimate the intensities of a particular 
characteristic X-ray line. The other is that, since peaks 
are about ten times broader than with diffraction methods, 
ten times as much of the underlying background intensity is 
incorporated within any one peak. Attempts to perform 
quantitative energy-dispersive analysis must include 
handling of these error sources as well as the usual 
corrections for absorption and other effects, as described 


below. 
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2.2 Quantitative Electron Microprobe Analysis 

As previously mentioned, the ratio of the intensities 
of an X-ray line in two samples is approximately the same as 
the ratios of their concentrations of the element concerned. 
More accurate analysis, however, requires the introduction 
of atomic number, absorption and fluorescence corrections. 
These are discussed in detail in Springer (1976). The atomic 
humber correction includes corrections for the back- 
scattering of reflected incident electrons, which increases 
With the average atomic number of the sample. Also included 
is the stopping-power correction which takes account of the 
decrease in electrons per unit mass as the atomic number 
increases. Thus unit mass of elements of low atomic number 
will absorb a larger proportion of the incident electrons 
than unit mass of elements of higher atomic number. The 
proportion of effective (X-ray producing) electrons 
increases with atomic number due to stopping power (S) but 
decreases due to backscatter (R) effects. 

After X-rays have been generated they may be absorbed 
by the sample before they leave the surface on their way to 
the detector. The ratio of transmitted to incident intensity 
is given by exp(-f,€t) where t is the distance travelled,€ 
is the element density and Bm is the mass absorption 
coefficient for that element for some X-ray energy Ey. A 
formula has been derived (Philbert, 1963), making some 
assumptions about the shape of the volume of sample excited 


by the electron beam, to express this ratio f£(X) as a 
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function of Ze As PryrOrE, and Eo- Z,A and Pm ate the atomic 
number, atomic weight and mass Ateop rude coefficient for 
the element, @ 1 Eo and Eo are the take-off angle, operating 
voltage and critical excitation energy for the analysis 
line. 

After generation, some X-rays absorbed by the specimen 
may have sufficient energy to excite the characteristic 
radiation cf another element present. The X-rays that excite 
the spurious characteristic radiation may have come from 
another characteristic line or from the background 
continuum. The corrections are different in each case and 
are rather elaborate. Springer (1976) summarizes current 
atomic number, absorption and fluorescence corrections, and 
Smith and Gold (1976) describe the equations chosen for 


correction of energy dispersive analysis data. 


3 Energy Dispersive Correction Procedures 

Because many corrections applied to X-ray intensity 
measurements assume a knowledge of the sample composition, 
corrections for samples of unknown composition must be 
applied in an iterative fashion, such that a previous 
estimate of the composition is used to apply the corrections 
to the raw intensities, therefore producing a new and, 
hopefully, better estimate of the composition. When 
successive estimates agree closely it is assumed that 


convergence on the correct composition has been achieved. No 


difficulties have been encountered with this approach if the 
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initial estimate was reasonable. For energy dispersive 
analysis in particular, correction procedures may be divided 
into instrumental calibration corrections, background 
shaping and scaling, peak overlap correction and composition 


estimation. 
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The first requirement for analysis is to know the 
energy associated with a particular X-ray line intensity. 
Although an instrument may be designed, for example, to have 
the lower edge of the first channel set at OeV and the upper 
edge at 10eV this will not usually be exact. The approach 
developed (Smith and Gold 1976) is to collect the spectrun 
of a calibration standard at the beginning and end of every 
analytical run. A calibration standard should give strong 
peaks towards the low and high ends of the spectrum. The 
common upper limit for analysis being 10 KeV, (10% electron 
volts), willemite was used as the standard, having low 
energy Zn Le and Si Ke peaks and high energy Zn Karp peaks. 

These calibrations permit the estimation of two sets of 
instrument parameters for the analytic run - those 
concerning peak position and those estimating instrument 
resolution. In both cases the first step is the 
identification of any peaks present and the precise location 
of their centres. Peak identification is performed by 
convolution of the spectrum with a 13-point filter as given 


in Savitzky and Golay (1964, Table III). The filter shown 
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gives the first derivative of a quadratic function (similar 
to the top of a Gaussian peak) and hence a change in sign of 
the derivative from positive to negative will give the peak 
location to the nearest channel. A cioser approximation may 
be made by making several estimates by interpolation of the 
full width half maximum (FWHM) of the peak and averaging the 
mid-points found. The result is usually correct to within 
0.1 channels. The FWHM, the usual way of expressing the 
width of characteristic peaks, is defined as the peak width, 
measured at a height of half the maximum peak intensity. 

The FWHM of a spectrum varies with energy (Smith, 


1976). The coefficient E, gives the FWHM at zero energy, and 


n 
the Fano factor F is a measure of the increase in FWHM with 
increasing energy. The FWHM at energy E is JE, 2+21*FE. Hence 
F and E, for the run may be determined from two peaks in the 
calibration standard and conversely, once E, and F are known 
the FWHM is Known for a peak at any energy. 

If the two energies of the calibration peaks are found, 
the peaks are identified and the true energies known, it is 
straightforward to calculate the instrumental 
stretch(scaling error) and shift (origin offset) for the run 
and correct all other spectra. The F and Ey, parameters are 
used to estimate the FWHM of each analyzed element, and 
hence the window over which integration of intensity is 
performed. In this fashion further processing of the spectra 
is independent of the instrument calibration. 


After a spectrum has been corrected for these 
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instrumental parameters, escape peaks are stripped from it. 
Escape peaks are produced within the detector itself because 
there is a smali probability that an incoming photon from 
the sample may excite the Si dead layer on the surface of 
the detector and that the resulting characteristic Si Ke 
radiation may re-emerge through the window of the detector 
and be lost. The correction for this effect is described in 
Reed and Ware (1972). 

Instrumental and escape peak corrections are applied 
only once to each spectrum. Other corrections, being 
dependent on the current estimate of the sample composition, 


are applied once for each iteration. 


Ware and Reed (1973) note that the detection limit in 
energy dispersive analysis is governed more by the 
limitations in background corrections than by the counting 
statistics. It is therefore of importance that the 
background correction should be as accurate as possible, but 
prior to the current work knowledge of the X-ray continuum 
shape was incomplete. Ware and Reed (1973) suggested that 
background intensity I at energy Ey may be calculated from 
the expression of Kramers (1923): 

[5.2] I = kZ(Eg-Ey) /Ey 
where k is a constant, Z is the average atomic number 
(proportional to the concentrations of the elements in the 


Sample) and Eg, is the operating voltage used for the 


Me 


Be Ss Toate ee 


* @ at a aleagie tS 
a mas pong ay by | Hee Liat ee 
suneried dime ate, job bat ee a7 eS ae 
oP Oe Mad oh — 7 
‘get? ans way eT) ne 53 shine Wo hom ithe “s 


Sete a2 orate 5 ba este vate 


pet? 4Pisy” Leto t eae ot Feluaba bai agi ti 
dt? to) Wel, oe Kothosi odgond- i 


soene sea St | 

= 7 j 

pew hs \ali Syeeee te eee eae nbs Beta e ac ee a 

= : m. ay 
. | / enti 

LI DSIVIOD. ALS Fy jeqstes hs besa 

Nsoa gkaROS 7 OG TT! sto, TI sth om aes i 


ane 290 sina 3 3 sea a 


> = 
a ae Te 


ey 
7 
i 
Ao 
al 
x 
ei 
, 


; ena. te be & fakes 302% a | 


aecoepga:l mdesseI79h 's29 7anD" 970i. 0G Very ae Nes ee 
att ; 7-7 ; 


som  Sokxstoo, ot e fal leas sv be ie 


i 


Sis) 9 


S42. va Nant wnetsoaT1e3) Sesospased i 2 


\ Pasi voo 


e 
‘ ay. 


: ”  gapu epee adaetiog ai \' 49) saoreie ads ae 


seh hha hae “Se *AT ODOe: 2h (9G ‘Debate doksoes oD. ha 


Lo 
ar 


cena sho verse ahr) 35 Se bat ides iol eee 


: bs a 


Jandy, L2reapple EXE boot bits aeK aotagane 


“2: bernie Lene rah se nee ie28k> SE Ll 
5 (ESE7 sniuban cor 


NEN IS gee = ‘a 
A aS | at weit otgote auageKs algae 2 ssn 


am! pont 42 a aieasis Rae x0 she Si 22; SE 


A ‘Silaraod bean aa i ines 


a 
aor 
a 7 eae 


cone 
q Ae { e i Ne a a. 
‘ f - ; ay : an 7 oon 
Rpm ve 


analysis. Additional corrections are made to this 
theoretical spectrum, primarily for absorption in both the 
sampie and detector. This spectrum is then fitted to the raw 
sample spectrum by scaling at two positions where peaks are 
unlikely. Correction for backscatter was not attempted, but 
Rao-Sahib and Wittry (+972) considered this to be necessary 
and supplied the relevant coefficients. They also modified 
Kramers" expression for backgrcund intensity to 
Poss ged,Ta= kyngseyoa? 

where x=1.11 and n is a function of E,, Ey and Z. However, 
as noted in Smith et al. (1975) I, as the number of photons, 
should be dimensionless, as is Kramers' expression, whereas 
Rao-Sahib and Wittry's expression includes units of energy. 
The precise definition of n was not given, aithough it was 
Said to vary between 1.2 and 1.37. 

Spectra were collected from the following substances: 
B(metal), Cq(diamond), MgO(synthetic), A1,0, (synthetic), 
Si0g (quartz), Si(metal), TiO, (synthetic), Fe (metal), 
Fe, 03 (hematite), Ni, Si(synthetic) and Cu(metal). Initially 
these were obtained at Eo = +5 kV kilovolts). Escape peaks 
were stripped from the spectra using the method of Reed and 
Ware (1972). Some of these samples were also run at 5, 10, 
12957520, °25 and) 29.5 kV) operating voltage. 

In order to use Kramers! background expression as 
Suggested by Ware and Reed 41973) it is necessary to correct 
oe absorption effects in the detector itself, especially in 


the low energy end of the spectrum. This requires accurate 
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knowledge of the thicknesses of the beryllium window, gold 
contact surface layer and silicon dead layer, with the 
consequent possibility of error. To avoid this problem Smith 
(personal communication, 1974) suggested that the ratio, 
(channel by channel), of two spectra of differing atomic 
number be used, with adjustments for any differences in 
probe current or counting time. Thus any deviation of this 
ratio from unity must be due either to peaks, where it would 
be ignored, an atomic number-dependent expression for the 
contimuum intensity, or differences due to incorrect 
absorption corrections. Various attempts were made to see if 
errors in calculating the atomic number and the absorption 
matrix corrections could explain the deviations of observed 
spectral ratios from unity. Mass absorption coefficients 
were calculated by the formulae of Frazer (1967), Heinrichs 
(1966) and Kelly (1966). The parameters h and fr used in 
estimating the actual percentage absorption of a particular 
Wavelength in a particular sample were varied widely. 
Backscatter corrections were made using the expressions of 
Duncumb and Reed (1968) and Rao-Sahib and Wittry (1972). 
None of these adjustments were capable of producing the 
desired unit ratio between pairs of spectra. Subsequent work 
was performed using Heinrich's (1966) mass absorption 
coefficients, Philibert's (1963) absorption corrections 
formulae and Rao-Sahib and Wittry's (1972) backscatter 
coefficients. Escape peaks were stripped by the method of 


Reed and Ware (1972). 
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Since the ratio variation could not be explained by 
errors in the matrix corrections, it became necessary to 
examine continuum intensity expressions. This was done by 
correcting each spectrum for matrix effects, as described 
above, dividing it by the proposed continuum expression and 
then dividing one spectrum by the other. If the sample- 
dependent component of the expression was adequate the ratio 
should be unity except for peak regions. For Kramers’ 
continuum expression - that is, using each spectrum divided 
by its atomic number - the ratios are approximately correct 
at zero continuum energy but deviate extensively with 
increasing energy and atomic number ratio. By using ratios 
of spectra obtained at one operating voltage the terms k and 
(E 9- Ey) /Ey in Kramers' expression, as well as any detector 
efficiency terms, are eliminated. 

Since Kramers' expression was inadequate and the 
correct expression clearly involved E,, Ey and Z2 an attempt 
was made to use Rao-Sahib and Wittry's expression involving 
Z . Since they stated that n was dependent on Z, Ep, and Ey 
but did not give the relationship, a simple additive model 
was attempted first, i.e. 

(5.4 ] n= at bEgj + cEy + da. 
Since ratios are calculated for spectra collected at the 
same operating voltage b cannot be determined from a single 
ratio plot, hence the expression is reduced to 

e5358 n= g + CEy + dZ 


where g = a + bEg. Examination of the ratio at Ey = 0 
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(obtained by extrapolation) showed this to be slightly 
greater than unity. A range of values of g and 4d were tried, 
and the best approximation was found with g = 1.06 and 

d = 0. The Z-dependent term in n was therefore dropped from 
further calculations. 

If spectra of two substances A and B, obtained at the 
same operating voltage are ccmpared, the two background 
Shapes can be expressed as Ae and n is the same for both 
spectra at any energy Ey, then the ratio of background 
intensities 

G54Gn a Tatpe= aelzzy! 
Hence it can be shown that 

564] n = log (Ig/Ip) /10g (Zg/Zp) - 
Fig. 5 shows plots of n for Cu and C spectra at 5, 10, 15 
and 25 kV operating voltage. It can be seen that, allowing 
for peaks and detector inefficiencies, n is linear. It 
should be emphasized that the only assumptions made at this 
Stage are that the continuum intensity may be expressed as a 
function of ase and n itself does not contain a Z term. 
Kramers" value of n= 1 plots along the base line, and the 
ranges of n found by Rao-Sahib and Wittry (1972) are shown 
as dashed lines. While n may fall outside these limits, it 
does not do so for their combinations of E) and Ey. Since 
the calculation of n is very sensitive to intensity 
fluctuations, care must be taken to obtain spectra at the 


same operating conditions and with as different atomic 


humbers as possible for the ratio Zq/Zp - In addition, 
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regions close to absorption edges should be ignored due to 
minor errors in the absorrftion corrections. 

The previously-described coefficient g may be obtained 
from Fig. 5. Its logarithm, plotted against the logarithm of 
Eos is shown in Fig. 6. Fig. 7 shows the slope of the lines 
in Fig. 5, plotted against 109, 9 (1/E5) - Both plots are 
approximately linear, and in this fashion errors in 
Obtaining these coefficients graphically may be reduced. 
When the adjusted lines are plotted on Fig. > it isrclear 
that they intercept at a common point (N,e in Fig. 5). This 
Suggests a further modification to the expression for hn, 
which may now be rewritten as 

(5.8 ] n= N + £ (EQ) (Ey/e-1), 
where £(E,) is the change in n for e keV and is a function 
of Ep. Fig. 8 shows a pLOt Of £(EO) against the logarithm of 


E showing good linearity with 


oO” 


eRe £(E = 0.2532-0.0584 1n(E 


o! 0° 
The resulting expression for n is 

[5.10] n = 1.159 + (0.2532-0.0584 in (E,)) (Ey/2.044 - 1) 
Which may be simplified to 

Poeitl) Wo =) tel 5o +) (0602 39=0./02857 In (E,)) (Ey~2-044). 

In summary, if characteristic peaks are ignored, a 

spectrum may be corrected so that its shape is independent 
of its composition, and hence of its average atomic number 
Z. This is achieved by dividing the observed counts by iM 


where n is defined above. If the shape is indeed independent 


of Z then the expression z" is a good estimate of the shape 


ot Sab er ones. ade DL gods note 


Lison re 9) So yay G $hb 25° T2409 h tae . 


yo) ehhbwayel weds se72505 i Te ae = aby 
22 euanees @irh “a Pee 


sO) GHP OLA AMS, Fe 


ALG es ni ; PORE ss cal \hy ores { cant 
tis ®% imme ie | enue ernie 
heaniees Gi tee Yhlsotdg ea 21 eae 

hie eee (688 He o Ste dahactan samt t 


Jon 
noir Oe eae! 


ue 


ry 


jc fet - tit a PAS re = wh) . (fug 


25 , dead ee el 
: vitt- “SAV (ene Fo Wisc A te 
ie _ a 

sobtomuatr att DAB, Von's Sum a as donndy coat 
10 ware yeeadi say Pantecs (oF ie ioigy s avoid & , 
ia piiaaes he eae 
» (ae me i200 sees, a = ages is 
ai a Lied adkoeouixe iy 

Bas Cea! OSES es on + eetigt a * ea 

ot sarretaniin het pa 


nds “4,31 weg Ten $0". Unt Sh Op ¥ deh wrt bine. 


a: isd 


7 
-_) 


s Story S3& ‘i isa) so | S249 S07 HAD, ay ae ou: Lies 

. - ay, 

7 atevnandouits al Saaht exh $psid OF, ‘sa sbaan “s La Sealed | 
= ie in 
ete pa oha9gts' sphas ve’ Pere a6 apiion bid ice . ase: 


ya qa etawor bowai@e do. ong (equine ‘gi, run 


apehancebat, baekns et ‘sgpda ie. ar son a 
Foal aa do etqmid ae pid” sane : mel ai 


n 


o 
ro) 


Ls Zo) He Ym NL 8 Ci Sel Seat Wa St fat Sa ak RY RATT A 2 Shek ont SY sy Sea La Se Hat at RT et TG} 


150 


140 


130 


—d 
NM 
So 


ro) 


3 3 | 

> H xe | x | 

YU H >| >| 
Og 


! 
fr ee 


CONTINUUM = ENERGY ( ke V ) ————® 


Fig. 5 (From Smith et ale, 1975) Plots of n versus 


continuum energy for various operating voltages, 
obtained from the ratios of copper and diamond 
spectra. 
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Fig. 6 (From Smith et al., 
relationship of f{Eo) to Eo. 


1975) 
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of the X-ray continuum. 

Fig. 9a shows scme of the ratios of spectra corrected 
using alle Agreement with the expected unit ratio is 
excellent, as opposed to Fig. 9b which shows similar piots 
corrected for n = 1 only, as suggested by Kramers. It should 
be noted that no fundamental significance is necessarily 
implicit in the form of the expression determined - other 
expressions may produce even better fits to the data. In 
addition, the number of significant digits in the 
coefficients for n may be fewer than given due to the 
unknown magnitude of errors in absorption and backscatter 
corrections and instrumental parameters such as probe 
current and operating voltage. The new expression does, 
however, provide a superior estimate of the continuum 
intensity than previous equations. This is indicated by the 
close approximation to the unit ratio in Fig. 9a, showing 
that those factors still affecting the continuum shape are 
not dependent on the compositional differences between 
Samples. 

In order to use the continuum expression directly to 
Simulate a background curve for a sample, it would be 
hecessary to know both the Eo dependent term -(EQ-Ey) in 
Rao-Sahib and Wittry (1972) - and the detector efficiency 
curve due to detector absorption of X-radiation. Hence a 
less error-prone method was implemented instead. Since the 
operating voltage of the system is not changed frequently, 


and the detector efficiency should remain constant, a 
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practical approach is to collect the spectrum of a standard, 
diamond for exampie, whose composition is known and which 
has no characteristic peaks in the region of interest. The 
spectrum is corrected for matrix effects and divided by vane 
This "normalized" background should therefore be free of any 
composition-dependent effects but should reflect any 
operating voltage and detector-efficiency characteristics, 
Once prepared, tha cauopuedeeed background is preserved and 
used for every specimen analyzed under the conditions at 
Which it was taken. Given the estimated or known sample 


a and 


composition, the background is multiplied by 2 
corrected for matrix effects. va is caiculated as (DiC; 20" 
Or, more accurately but with greater computational effort, 
as Licj i) where Cj is the concentration of component i. 


The resulting spectrum should have the shape of a background 


produced by a sample of that composition. 


—— = setae oes eo ew ee ae ee Se we oe 


Since different spectra will be acquired with different 
probe currents and counting times, the shaped background 
must be scaled to the sample spectrum. Although this could 
be done using the probe current and counting time, it is 
more accurate to base the calculations on observed 
intensities, as this eliminates error due to instrumental 
drift. Ware and Reed (1973) suggested that scaling be 
performed using the Ar Ke and Ge Ke energies (2.96 keV and 


9.876 keV respectively), which will not normally be 
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interfered with by other peaks. Disadvantages with this 
approach are that the Ar Kqe position may be interfered with 
by residual pulse pile-up from Al and the intensity of the 
Ge Ke position is very low, producing poor counting 
Statistics; 

In order to scale the whole spectrum to an overall best 
fit, it is first necessary to eliminate peaks in any 
averaging. As a first step Rye the lowest spectrum over 
background ratio (performed uSing a 13-channel running 
average to reduce statistical fluctuations) is obtained. If 
any of the background is clear of peaks, this ratio will 
approximate the correct scale factor for the shaped 
background. An upper limit Ry is then calculated for each 
channel in turn. A final average ratio is calculated using 
all channels whose intensity (cr number of pulses) I is less 


than R thus eliminating all significant peaks. The upper 


u’ 
limit is calculated as 

[5002 Japk yrs Ry (143. 5/,/T) « 
The reciprocal of JE is the standard deviation of the number 
of pulses I, if I is sufficiently large that the binominal 
distribution of the counts approaches the normal 
distribution. Thus the upper limit of the acceptable 
background ratios is set at 3.5 standard deviations above 
the lowest ratio found in the spectrum. The number 3.5 was 
found empirically to provide the best discrimination between 
peaks and background for the instrumentation used. Energies 


below 2 keV or above 9 keV were not used, due to the steep 
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slope of the continuum at low energies or the very low count 
rate at high energies. All spectra were adjusted for stretch 
and shift, as described earlier in section 5.4, prior to 
background fitting. Fig. 10 shows the effect of background 
subtraction on the spectra of ten typical substances. In 
general the procedure is extremely satisfactory. Slight lows 
just to the high energy side of some peaks are probably due 
to inadequacies of absorption corrections for characteristic 
radiation when applied to the continuum (Smith et al., 1976). 
Residual counts on the low energy side of peaks just above 
Si in energy, are iow energy tailing and shelving effects 
due to incomplete charge collection in the silicon detector 


(Smith 1976). 


Since energy resolution is much lower in energy 
dispersive analysis than in wavelength-dispersive analysis, 
the overlap of peaks is a frequent occurrence. While 
attempts have been made to distinguish peaks by assuming 
Gaussian shape, many overlap conditions are not amenable to 
this approach, particularly where the non-Gaussian tails of 
peaks are involved. Consequently a more general approach was 
applied. 

Peaks are separated so that a central part of each may 
be integrated to produce a measure of the total intensity of 
that K-ray line. Since F«and E, are estimated for the run, 
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and spectra are therefore automatically corrected for 
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stretch and shift, it is possible to define a run-condition 
dependent set of integrating regions each covering 1.2 fuil 
width half maxima around the peak centre for each analytic 
line. It is therefore possibie to process a set of simple 
elements or compounds and obtain the above-background 
intensities that fall within the defined integrating 
regions. These are expressed as a fraction of the intensity 
obtained for the analytic line of the element being 
processed. 

Table 5 shows the overlap coefficients for 22 elements 
commonly used in analysis of silicate minerals. Fig. 11 
Shows the main features of this table diagrammatically. 
Overlap coefficients decrease progressively on the low- 
energy side of parent peaks, but increase as the energy of 
the parent peak approaches that of the silicon absorption 
edge. "This pattern indicates that the shelving effects 
reflected in these overlap coefficients, result largely from 
partial loss of charge in the silicon deadlayer, as well as 
from trapping in impurity centres in the detector." (Smith 
and Gold, 1976, p. 196). In addition, where overlap is 
produced by element peaks other than the analytic line, 
differential absorption corrections should be applied. This 
is the case for the Ka peaks on the high energy side of the 
analytic Ke line. In addition, atomic number and 
fluorescence corrections are applicable where Le peaks 
overlap some regions cf interest and analytic results are 


obtained using the Kypeak for that element. 


ra 


Di) I<6 $0, Yee s) 2eeodee,. OF anaes 
pidoaptined-— Woe) ds aLasio bas et its 

par psty aris bu etd GF eo Li 

yoienetat Sit to sohaa: 3 i 


salted. thane ls ads 


he 
—, 
ad 
6 —= 
j 
va 
i 
;> 
> 
w 
rm 
a 
t 
r 
— 
hw 
4 
rw 
~~ 
4 


fy? L648 wabararin stsotl fe to 


<vipd od 10) YASV EGE S202 $5 baIseo 4 

10 RO hier. ou eaeaq) 3 
sd20eGe, nogt lta alt 25 rhie, aansanniae aia 8 

pros tievpsivteds civ rodt escepebay vantaon kde . 
aet YlepreL Tinses «2 lisioiatags) oa niki a 
a ira imes fa VBE fob ndvkiia Sat ‘kt sail seman 
atime) " OED Seb ode one anranas «ot aAiyil ae 2 was 


ef 1aiuev SISA. | ANRGLSD a at (2ST at aie 
elid, srsyleus, pas tad ¢. isito-shead ana. 


tp 


gan! ,f6ecis92 54 pdudiie, BiObFos 109 age "a 


im 


Su? 20 ebie yorsns det ad +" ap rice Ae aad ate 


- Bids isagta bei “\dobvithe r: pe rin 
akeqa~ gi 9292M 2 sdnoniys 225 “socgabaa scape goa ni 

S16 271GAS3 sires 6a8 bite, Fe aweg dt! ia ane 
‘asfels Fad zo¥ Yoeley ae 


—< 


a3 


“UOLYNLOSAA 4OQIaZaP ueyNdt 
quedijiubis Bursnpoud aqjeu Burqunod uw: sabu 
3e AluO *A[YDLUzs Spue (Ay GL) pesn abeqLo 
SUOLJIOUNOD {YZ LeLJUSUSJJLP L [NZ ‘SaULL y PazeLoOSse jo sabequacuad se Ppessaudxa aue squ 
"S2uaueaja Gurdde(4aao yo syead sisfjeue ut (WHM4S Z2°L) S®tatsuaqut paunseaw yo sabequas 


ee 


O01 12k 
; Ore mns 
001 40°0 20°0 
1Z°0 OOL OL'O 90°0 €0° 
Of"L OOL EL°0 60° 
* 0S°2 OOL SL'0 
e 00°S 00 
e 02° 
= 21°0 
Eat 
= OL'€8 
a 
o 
a 
S Li'l 
2 pb'e 
1°88 


qued e@ 03 ALuO ajqediidde ase Aauy *OSLy 


eu) 


6L°2 66°2 


S9°el 


(9L6L 


BEEE €9°S6 £0°2 


c0°0 

90:0: 20-0 

€t°0 90°0 

92°0 SLO 

OOL Z2é‘0 

OZ°€L OOL 
SE-cl 


€O°OL 00 
L° 


quauea,[s paddeiusaao 


Ae 


c0°0 
£0°0 £€0° 
£L°0 
0v'0 


‘pTop pue yytms word) G¢ eTqel 


SLU2 UL S3uaLdis4a09 


£50" £0<0 


oot 


MYM-—-OO 
Cooooo 


@ 
° 


WwW Ne © 
cCoooeo 
TMOM WO — 


wo 

nC 

=) 

pa pe 
-—-OFSNnN~O 
oOn—-- oO 
CoN (coy Cas (cm iC} 


MO DWr- ON 


im +t iwoOor- CO © 
Seer CS CCC 


fon) 
—- o> 
a) 
ore 
NS 
- © 
S| 


°S$ZUSLILYJBO9 |yZ ohueyd [LLM HuLuapeoug yead 
*Peulezqgo suam Aay 2 YydiumM ze (unuzoeds LINZ) Das/Q00E JO 2224 BHulquNoOd abeuaae ayy 
A Bulqesado ay} ye Aluo ayqeoijdde ase ajqez 


“atqediidde oue 
7"O°N 


PLILLJ9OD de_Asao sul, 7 auayM 
dad se passaudxa sjuatdiyya0d deysaag 


80°L 
69°0 
92°R 
BLL 
bel 
Sea 
62°71 
cc°0 
S2°0 
o¢ °0 
ve "0 
8t'0 
Sv'0 
60°0 
S520 
02°0 
96°0 
62°0 
Se°0 
ene 8) 
LO'2 
00l 


u| 


quawa,ea BuiddepusrAo = ----- 


a 
as 


& Weil 


“Es? > 


a8 at 


1 
tel we a 
» =, 1 Pasi 
7 ati + * A , = 7 ao 
| ae a: 
™ : 4 “Ff ; 7 2 ; 
= »* ; 
i. 9 


lying) pol race ty 
od ae 
Mane 


_ 


a) ee 
ae 
ee 


54 


“G eTqeL jo uz 
04 ad B8O1IZ Sqyuatotryzsoo detisao ayy Jo uoTtpezueserdsarz 
OTJeuBeriberq (9/61 “PpTOD pue Y3THS woz} iis -*Dta 


Lipp dVIIAO INVDISINOIS yy 
Litt ee ON JO NOID3¥ Gy) 
NU ean.” oe 

aa 2 


dVIAFAO INVOISINOIS 
ON JO NOI9S4 


—1IN3W313 ONIddV183A0 —> 


& 
~N 


> 
. 


"au 


a oe 


> 


re i 


7 
iu 


% 
Tt 
i a 


e 


Since the measured intensities for a sample are the end 
result of all relevant overlaps, and the overlap 
coefficients are known, the solution of a set of 22 
Simultaneous equations for the 22 intensities measured at 
the locations of the analytic peaks provides the solution 
for the 22 original intensities prior to overlap. Overlap 
coefficients are, however, affected by the system 
resolution, which is affected by the counting rate. Hence 


all spectra should be obtained at similar counting rates. 


2:8 Calculation of Composition 

Once intensities are obtained they are corrected for 
matrix effects as described in section 5.2. Corrected 
intensities obtained from a sample of unknown composition 
must then be compared with intensities obtained from 
standards of known composition in order to estimate the 
concentration cf any particular element. Thus 

5.135) Cgp/Isp = Cst /Ist 

where sp and st refer to specimen and standard, 
respectively, and C and I are concentration and corrected 
intensity, respectively. Standards must therefore be 
processed in each run for each element analyzed, prior to 
the analysis of an unknown. Fig. 12 shows a typical printout 
for a standard, using the EDATA program (Gold and Smith, 
1976). While this is undoubtedly the most accurate method, a 
considerable saving in analysis time can be made if 


approximate standard intensities can be preserved from 
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Fig. 12 Sample EDATA printout for a standard. 
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previous runs for use in estimating the concentrations of 
minor or unexpected elements. For this purpose standard 
intensities are preserved in a standardized form suitable 
for use in any run. Intensities are recorded as counts per 
second per nanoamp of probe current per percent of that 
element present in the standard. While not as accurate as 
run-time standards they have appreciable value. Fig. 13 
shows the EDATA printout for a typical sample analysis. In 
the last line under "OLD STDS.", can be seen the standard 
numbers used for each element. In this list a zero indicates 


that default standard values were used. 


5.9 Seguence and Implementation of Correction Procedures 
The normal procedure for analysis of an unknown 
therefore contains the following steps. 
1/ 
The sample spectrum is corrected for stretch and 
shift and escape peaks are stripped from it. 
2/ 
The spectrum is scaled to standard probe current. 
3/ 
An initial estimate of composition is made. 
7 
The normalized background is shaped using the current 


composition estimate, involving the application of 


2 n : 
Matrix and Z corrections. 


5/ 
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13 Sample EDATA printout for a sample of unknown 
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The background spectrum is scaled to fit the sample 
spectrum, and subtracted from it. 

6/ 

Integrated intensities are obtained over standardized 
regions of interest for each element, using the 
background-stripped sample spectrum. 

is 

The integrated intensities are corrected for peak 
Overlap. 

8/ 

Resulting intensities are corrected for atomic 
number, absorption and fluorescence effects. 
9/ 
New compositional estimates are made by comparison 
with standard intensities. 

10/ 
If the current and previous composition estimates are 
sufficiently close (within 0.05% for all elements) 
calculation ceases. Otherwise the procedure is 
repeated from step 4. 

Fig. 14 shows a schematic flow chart of the EDATA 
package. After subroutine INIT initializes parameters for 
the run a spectrum is read and standardized by subroutine 
READIN, and, if this spectrum is of a sample of unknown 
composition, the initial composition is estimated from the 
raw peak intensities by calling SCAN with a parameter of 


zero. After this , subroutine SHAPE is called to shape the 
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14 Schematic flowchart for the EDATA program. 
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normalized background to the sample spectrum, SCALE with a 
parameter of zero is used to fit the shaped background to 
the raw spectrum, and SCAN with a parameter of one is used 
to search for and display all major peaks and their size 
with respect to the background (see Figs. 12 and 13). 
Following this, SHIFT is called to adjust the raw spectrum 
for any instrumental stretch and shift as estimated from the 
calibration standards, ESCAPE is called to strip off the 
escape peaks, and SCALE with a parameter of one does a final 
scaling adjustment to fit the background to the shifted 
Sample spectrum. Subroutine INT then integrates both the 
shifted spectrum and the scaled background over the region 
of interest for each of the 22 elements, as estimated from 
the FWHM of the calibration standards. The spectrum-minus- 
background figure for each element forms the intensity value 
used to generate subsequent compositional estimates. 
Subroutine TEST handles all subsequent processing of 
intensities. OVLAP is called to estimate the original 
intensities in the absence of any peak-to-peak overlapping. 
ZAF corrects these intensities for atomic number, absorption 
and fluorescence effects. Processing of standards terminates 
at this point with the derivation of standard correction 
factors in the form of weight fraction of an element per 
ZAF-corrected count / second intensity, to be used for 
compositional estimates on subsequent samples. Thus, in Fig. 
12, the values in the "STD." column are obtained by dividing 


the known standard composition values in the "WT.FR." coiumn 
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by the intensities in the "ZAF*C/S" column. In the case of a 
sample being analyzed (Fig. 13) the new concentrations 
q"wl.FR." column) are calculated by multiplying the 
intensities in the "ZAF*C/S" column by the values under 
"STD.". Subsequently a user-defined AJST routine is called 
to modify the concentration estimates, e.g. by estimating 
the oxygen concentration as the difference between the 
compositional total and 100%. This compositional estimate is 
then used in loop 2. On the second or subseguent iteration 
only SHAPE (to reshape the background on the basis of the 
new composition estimate), SCALE with a parameter of one (to 
scale this background to the initial shift, stretch and 
escape peak corrected spectrum) and INT (to calculate peak 
intensities) are called prior to re-entering the TEST 
subroutine. After new concentration estimates are obtained 
they are compared with the previous values, and a further 
iteration is performed if any element deviates from the 
previous value by 0.05% or more. If, however, convergence 
has been achieved, the user-defined subroutine CALC is 
entered (perhaps to recalculate elemental concentrations as 
oxides), SCAN with a parameter of two searches for peaks on 
the background-subtracted spectrum (Fig. 13, bottom) and the 
program proceeds to read and process any subsequent spectra. 
Further details of the computer implementation of this 
approach as the EDATA package are contained in Gold and 
Smith (1976). Typical analyses obtained using this system 


are shown in Table 6, from Smith and Gold (1976). For major 
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elements the accuracy is possibly better than that 
obtainable by wavelength dispersive analysis, since all 
elements were analyzed at exactly the same points on the 
sample, and their intensities were collected simultaneously. 
For minor elements, present in concentrations greater than 
about 0.25%, accuracies are similar, but below this 
accuracies become appreciably worse than with wavelength 
dispersive analysis, although with care concentrations of 
500 ppm may be detected. Currently the detection limit is 
considered to be 300ppm. In Table 7 the best available 
analyses found in the literature (column 2) for an obsidian 
sample are compared with an analysis using energy dispersive 
analysis and EDATA (column 1). Elements underlined were 
calculated using standards processed during the run, while 


Ci and Mn were calculated using default standards. 


wa) 3 . : } mt - a Bis. ‘ : oy : 
sear meds sino iteseq at’ 
ia ae man sn 
es ones. enous ove ra pit i 


a 
Ce 


ode xo a sit tog ae B24 4h 


wvipwoe net iim. 2 Fs haateie ean sina 


aaby retsemy acokdss ine aaoNpa se: 


eae? aqued eUdl< 1h hee re oie. 


ae wae bes ad t,t De apa ay x Sit 


49 Byaeeort ke oo) Ss HS aE ha cnapa 


eo 


at Piapsk soteoessh sft ¢ eorasaup {oat fi ) 


eitelgave cea) ant \ akaet nt rr ba tare 
49 cn a 


. 7 Pia 


atinseqatb ypyeus. nofen.2isgians as aati 


sek biaide ap tot (es Hebloo) saps 
eet bong isiou sonemals, . (t medooy | Qa Dict, Se 

ja, nF ae 7 
. See 


etki eaua bat jasa0B aS shipbrede lay Pf 
Cp ifeekeh pitas ki per 


7 
7 


eS 7 ne Pe any ay m 
a | 7 


Table 7 Comparison of analyses performed by the 


progran 


time standards are underlined. Oxygen is 
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Contour Mapping Techniques 

Contour mapping is an attractive approach for the 
display of the chemical variation in the analyzed till 
samples, or the elevation of the contacts between 
stratigraphic units. Available computer techniques were 
considered unsatisfactory for a variety of reasons. Computer 
programs using traditional surface estimation methods to 
produce contour maps follow a fairly Standard procedure. A 
Suitable rectangular grid is superimposed over the map area, 
and an elevation estimate made at each grid node. Estimation 
techniques vary widely, but usuaily consist of selecting an 
neighbouring data points, evaluating an inverse distance 
weighting function for each (giving less weight to more 
distant points), and then calculating a weighted average for 
the values obtained from these data points (€.g. Sampson, 
1975). Once the grid node values are obtained, straight line 
(or, occasionally, curved) contour segments are obtained for 
a desired contour value by interpolation from the values at 
the corners of each rectangle. Reviews of the multitude of 
contouring methods are given by Rhind (1975), Walters (1969) 
and Crain (1970). Some of the better known publications on 
the subject are Shepard (1968), Cottafava and LeMoli (1959) 
and Dayhoff (1963). 

Several problems remain with most gridding methods. One 


is the use of an assumed arbitrary distribution function 
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such as 1/d#, where d is the distance of the data point 
from the grid point whose value is to be determined. An 
excellent discussion of distance weighting functions is 
given in McLain (1974). Another difficulty occurs in the 
definition and collection of the "neighbouring™ data points. 
This does not necessarily mean "closest" since it is 
desirable that the points surround the location under 
examination. Coliection and testing of the data points is 
frequently expensive since there is no direct method of 
locating the nearest points to any arbitrary map location, 
and so some form of search of the whole data set is 
necessary. If the weighting of a data point has not reached 
zero when it is rejected from the set of neighbours being 
used to estimate elevations in a region of the map, then 
discontinuities in the surface may occur with weighted 
average techniques. Further problems are (1) the grid size 
defines the final resolution irrespective of data point 
distribution since further processing consists of 
interpolation between grid nodes, not data points, (2) 
interpolation of contour segments in a rectangle is not 
always unique or trivial due to the difficulty in 
distinguishing a north-east trending ridge from a north-west 
trending valley (or vice-versa), (3) fairly extensive 
progranming may be required in order to string ali the 
contour segments together for plotting and (4) elevations of 
the defined surface may not pass precisely through the data 


points. For topographic maps this last property is 
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undesirable, but for data with an appreciable error 
component some smoothing may be desirable. Gold (1977) 
discusses some aspects of the selection of a suitable 


contouring method. 


Contour mapping using triangular element data 
structures possesses several advantages in dealing with the 
above problems. The general approach is to break the map 
region into triangles with data points at the vertices, 
obtain the slopes at each vertex and bend triangular plates 
to conform to these. Each of these irregular triangles is 
represented by an element in the data structure, with 
pointers to the three adjacent triangles and the three 
associated data points (Fig. 15). A data structure of this 
type that was defined manually rather than automatically is 
described by Heap (1972). 

In order to understand the benefits of the suggested 
data ee eaerure it is important to understand its 
construction and properties. The fundamental concept is the 
availability of a local homogeneous coordinate system rj, (1 
= 1 +o 3) associated with each triangle (Fig. 16a). Thus an 
arbitrary point p at location (Xn 0 Yp) may be located with 
reference to a particular triangle with vertices at the 
POIntshar=o1,57 2,. 3 (numbered anticlockwise) by 
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Fig.15 
structures. 


Data elements 


Triangle # 


Adjacent 
Triangles 


5,—,4 
3,6,4 
7 ea 
7 he 
1 ,6,— 
2,—,5. 


Associated 
Data Points 


1,2,5 
2,6,4 
3,4,6 
1,6,2 
7, 5,2 
72a 


(elevation ) 


(From Gold et ale, 1977) Triangular element data 
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a) Local coordinates 


i 
. ry 4 0-4 
9 = 0.2 
fy = 0.4 


& =Q-2 emesccmner 


Beer 13; v1 


b) Planar interpolation of elevations 


PA (P)= r & (Vy Re r9 Z(V> )+r3Z (V3) 


c) Local coordinates and search strategies 


Nr2 2-1 . 
r250 oR Point ory tg 8g 
~ vy han Rees Fel B LOA. 0.40.2 
P2 08 -0:3 0-5 
P3 -0-4 -0-6 2-0 


eo ry =O 
a i 
Seen ry s-]} 
f ? / 
/ / 
vd / / 
7 if 
r320 Lc =] r3 22 


Fig.16 (From Gold et al., 1977) Derivation and use of a 
local homogeneous coordinate system for a triangle. 
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The value dj may be thought of as twice the area of sub- 
triangle i, and D as twice the area of the reference 
triangle. This is the solution of the system of linear 
equations: 

(6.3] X= TY, X1 tly Xp try x, 
Une Gi" Peet a os 


e=tr +L 4 


; 2 3 


In addition, if planar interpolation of elevations is 
required, 

(6.4 ] Zp= Y Z, tlp Zp +1, Z, 
as illustrated in Fig. 16b. 

This coordinate system was first devised by Feuerbach 
(1827) for tetrahedra and variations on it were extensively 
Used by geometers of ‘that century. Each of these three 
coordinates (Ij ) has the value 0 when p falls on the base 
dine i, and 1 when p falls on a line parallel to base line 
1, but passing through the opposite vertex i. The three 
coordinates sum to unity, hence any two define the third. 
All three coordinates are positive if p falls within the 
triangle, and if p falls outside, the signs of the 


coordinates indicate the direction to go looking fora 


circumscribing triangle. 
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dieie clear from Fig. 16c that, if a data structure 
relating adjacent triangles exists, a search algorithm is 
easily devised to start at any arbitrary triangle and search 
through a set of adjacent triangles until the one enclosing 
the unknown point is located. This in itself suggests a 
method for generating such a data structure, A “universal" 
triangle, iarge enough to enclose the map area, is first 
generated. The data points are inserted one at a time and 
the enclosing triangle for each is located. This triangle is 
then divided into three new triangles with the new point at 
a vertex common to all three, and all pointers are then 
updated (Fig. 17a). 

The resulting mesh of triangles is continuous over the 
Map area but the triangles produced may have undesirable 
Shapes. An optimization routine is then used to test each 
triangle against each of its neighbours. The quadrilateral 
thus formed, if convex, may be divided in two ways, and it 
is the purpose of the optimization criterion to choose one 
of them. If the current division is the poorer of the two 
then the quadrilateral is divided in the aiternative manner 
and the data structure pointers updated as in Fig. 1/7b. In 
this manner every triangle pair in turn is examined to see 
if it is acceptable and, if not, the alternative partition 
is taken. This procedure is repeated until an entire pass 
through the data structure produces no changes. By the 


judicious setting of flags the optimization of even a few 
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a) Inserting new data points 


Map elements Data elements _ 
Before 7 Triangle # Adjacent 
Triangles 
9 1, 3,8 
9 1, 10,11 
10 9,3,11 
H 9,10, 8 
P4 
b) Quadrilateral optimization 
Before Re 
I 7 9, 5,8 
9 7,6,2 
Pg 
7 7 Mads ye) 
9 8,6,7 


P4 


Fig.17 (From Gold et als, 


Associated 
Data Points 


7,2,4 


9,2,4 
7,9,4 
7,2,11 


2 ,8,5 
4,5,8 


2,4,5 
4,2,8 


1977) Triangular data 
structure generation and optimization. 


16 


= 


* ss 

arb tM 
ia 
ee) 


whee | iver 
238% 


ie 
i ; - 


_ pat 
a ier 


74 


thousand triangles remains an economic proposition, the 
insertion and optimization of 1129 data points requiring 
0.22 and 1.92 seconds respectively on an Amdahl 470/V6 
computer for a test data set. These flags are set to define 
triangle pairs that have been tested and found satisfactory, 
and hence need not be tested again until one of the two 
triangles has been redefined. Fig. 18 «solid line) 
illustrates the rate of convergence of the iterative 
Optimization. The natural logarithm of the number of 
Switches made decreases by close to one third on each 
iteration, and hence the number of iterations approximates 3 
in(t) where T is the number of triangles. This is a worst- 
case condition since the dashed line shows the convergence 
rate for a second batch of data points inserted after the 
first set had been optimized. This approach seems to have 
considerable economic advantages over the methods of Frenkel 
and Gill (1975) and McLain (1976) in which the mesh is built 
up agglomeratively from some nucleus by testing all 
neighbouring points to find the one most suitable to be 
added at a specific mesh location. Lawson (1977) generates a 
mesh by first finding the convex hull of the set of data 
points and agglomeratively building the mesh within this. 
Optimal triangulation (as defined in Section 6.4) is 
achieved by testing triangle pairs, as is done in this work. 
The use of a fast triangulation algorithm, as suggested 
here, may permit the contouring of moderate to large data 


sets several times faster than by gridding methods. 
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Fig.18 (From Gold et al., 1977) Rate of convergence 


optimization of a triangular mesh. 
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An obvious question remains - what is our optimization 
criterion? Rhind (1975, p, 297), in reviewing spatial 
interpolation techniques, states his ignorance of any 
",..fundamentally sound, cheap and easily-applied 
partitioning theory". He does however, mention an Israeli 
firm's approach of a single pass through the lattice to 
Minimize the total length of the triangle sides, as well as 
a commercial package (SCA, 1970) which use the ".,..best 
Overall approximation to equilateral triangles" (Rhind, 
1975, pe 297). Gold (1976a) similarly optimizes by choosing 
the partition of the previously mentioned quadrilateral with 
the shortest diagonal. This however proves to be somewhat 
unsatisfactory as triangles with an unnecessarily small 
height could sometimes be produced. A more satisfactory 
criterion turned out to be the maximization of the minimum 
triangle height. This height is measured perpendicular to 
the diagonal bisecting the quadrilateral. Thus two heights 
are obtained for a given partition of the quadrilateral, and 
the smallest of these is retained. This height is then 
compared with the equivalent figure for the alternative 
partition of the quadrilateral, and the final partition 
Chosen is the one having the largest value for this minimum 
height. Fig. 19a illustrates the difference between this and 
the shortest diagonal criterion. The minimum diagonal 
criterion produces the left-hand triangle: however, the 


height of this triangle measured perpendicular to the 
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a) Differences between optimization criteria 


minimum diagonal partition 


maximized minimum height partition 


minimum height 


b) Imperfect optimization 


Fig.19 (From Gold et al.- 1977) Optimization criteria. 


(a) shows the difference between two proposed 
Criteria, (b) indicates areas for which triangulation 
is not optimum and (c) shows the geometric 


relationships between two proposed criteria. 
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relevant diagonal is the smallest of the four possible 
triangles, and hence this partition is rejected when the 
maximized minimum triangle height criterion is used. Lawson 
(1977) maximizes the minimum angle, instead of the height as 
is done in this work. 

McLain (1976, p. 179), referring to Pitteway (1973) 
states that "An optimai partition is one in which, for any 
point within any triangle, that point lies at least as close 
to one of the vertices of that triangle as to any other data 
point." He then gives an algorithm to choose from among the 
candidate points. While the definition is excellent it 
cannot always be achieved, as shown in Fig. 19b. Whichever 
way the guadrilateral is partitioned, one of the shaded 
zones will not be “optimali" since any point within that zone 
will be closer to a data point outside the triangle in which 
it occurs than it will be to any data point forming a vertex 
of the containing triangle. Thus any point in a shaded zone 
is closer to a "bad", or undesirable, data point than it is 
to any "good" data point. The worst possible situation 
eeecursvdtepotits Peor' Oi an Fig.u19cpvatithe closest-potntein 
the shaded zone to the “bad" data point. 

With this situation in mind Pitteway's optimal 
partition may be rephrased as an optimization criterion. 
This is defined as minimizing the ratio of the distance to 
the closest "good" point over the distance to the closest 
"bad" point, measured from the worst-case position in the 


triangle. If the closest "good" points for P and Q are C and 
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B respectively then the relevant ratios are PC/PB and QB/QC, 
or cotQ@ and cot@ respectively. Since simQ= PB/BC and sin@ 
= QC/BC it is clear that maximizing the minimum height (i.e. 
choosing the larger of PB or QC) fulfills the theoreticai 
optimization criterion. If the nearest "good" point is not 
on the suggested diagonal then PC and QB form upper limits 
to PD and QA, again fulfilling the required optimization 
criterion (Gold, 1977). 

It is, of course, sometimes advantageous to modify or 
build the triangulation by hand. This provides the facility 
to modify the surface on the basis of additional 
information, such as continuity of valleys and ridges, which 
is normally difficult to build into the mapping program. A 
good description of this problem under the heading of 


"Surface specific points" is contained in Peucker (1972). 


Once the data structure has been generated, its 
benefits in contour mapping, for example, become apparent. 
The problem of identifying the neighbouring data points for 
any arbitrary location is readily solved. Its three 
immediate neighbours form the vertices of the circumscribing 
triangle and, if desired, the three additional vertices from 
the three adjacent triangles may be included. This approach 
forms a very happy combination of the requirements for 
neighbours to be both close to and surrounding the point 


under consideration. If a simple contouring method is 
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sufficient, each triangle may be considered to generate a 
function consisting of a plane passing through the three 
data points, as in Bengtsson and Nordbeck (1964). The 
surface made up of the set of triangular planes is readily 
contoured, there is no difficulty defining the domain in 
which a given set of data points define the surface, and the 
resolution of the final map is a function of the local data 
point density. The surface is, of course, discontinuous in 
the first derivative across triangle boundaries and its 
contour lines are therefore angular. For many applications, 
especially where the data point distribution is very dense 
and/or the final product required is the elevations at the 
nodes of a grid, this is an extremely efficient technique. 
Alternatively, conventional contouring techniques could 
be used to produce a grid of values, merely using the 
triangular element data structure for collection of 
neighbouring data points. This would, however, neglect 
several advantages of the triangle approach. Therefore, in 
this work a method has been devised to interpolate within a 
triangle using the local coordinate system of Fig. 16a. In 
this approach, the first step is to estimate a best-fit 
plane, (constrained to pass through the data point), at each 
datum location. The surface is then estimated for each 
triangular domain in turn using only the three planes at the 
vertices. These three planes may be expressed ina form 
compatible with the triangular coordinate system as Bij , the 


elevations of each of the three planes i at each of the 
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three vertices j (Gold et al., 1977). If field measurement 
of the slope at data points is available (for example with 
dipmeter readings from drill holes) this may be used instead 
of the estimated best-fit plane. Use of Slope information 
greatly increases the information content of the surface. It 
is suggested in McLain (1974) that the use of x and y slopes 
as well as the elevation information at n data points is 
approximately equivalent to 2.5n data points with elevations 
alone. This approach was also used in the SURFACE II 


contouring package (Sampson, 1975), although in both cases 


interpolation was made ontc a square grid. 


Within each triangle an equation can be developed using 
the local coordinate system and the nine values Bij fora 
Surface possessing certain specified properties. These 
properties are similar to the surface properties for any 
contouring method. The first is that the surface passes 
through the three data points - a condition not met by many 
mOving average and gridding methods. The surface should also 
possess the measured or estimated slope at each data point. 
This condition is rarely met in conventional contouring 
programs, and with moving average techniques the slope at 
data points may be zero. It follows from this that extrema 
Should be free to occur between data locations. This is the 
case with global polynomial methods, although the values may 


be unreasonable. The SURFACE II program, (Sampson, 1975) 
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achieves flexibility of extrema locations by weighting the 
estimate produced by a best-fit plane at each data point. A 
fourth requirement is that the surface is continuous between 
domains, Or regions of the map in which a specific function 
is used to interpclate the surface. For triangulation 
techniques this domain is obviously defined as the triangle. 
For conventional methods the definition is less obvious, but 
in fact it is the region of the map where a particular set 
of "neighbouring" data points are used to estimate the 
surface. Failure to give a point a zero weight before it 
leaves this set results in unexpected breaks in the surface 
and kinks in the contour lines. A last requirement is that 
the slope of the surface is continuous between domains, 


providing smooth contour lines across domain boundaries. 


Several attempts have been made to meet ali these 
criteria, especially those of continuity of slope. 
Suggestions for obtaining an interpolated elevation estimate 
by weighting each of the three planes at the vertices of the 
triangle by the square of its associated local coordinate 
{C-value) (Gold, 1976b) or weighting a quadratic curve at 
each vertex by r3 (McLain, 1976) have been made. Both 
claimed that the slope of the interpolated surface did not 
Change abruptly when triangle boundaries were crossed, and 
in both cases this appears to be incorrect for reasons given 


by Zienkiewicz (1971, p.175). Differentiation of either 
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Gold's or McLain‘ts equations in Cartesian coordinates, 
although tedious, clearly shows the misconception. For 
example, the function slope at edge 1 of a triangle is 
indeed independent of the coefficients of the plane or 
quadratic surface associated with the opposing vertex 1. 
However, the function slope is not independent of the 
location of this vertex since the weighting functions used 
On the remaining two surfaces associated with vertices 2 and 
3 are dependent on the locai coordinates ly and tr which 
have been shown previously (equations {6.1] and {6.2]) to 
involve the values x and Y,° If moving vertex 1 changes the 
slope at and perpendicular to side 1 then clearly slope 
continuity cannot be expected when passing into the adjacent 
triangle. The effect is particularly marked when angle 2 or 
angle 3 of the triangle is obtuse. If a function (in this 
case surface slope) is to be continuous across triangle 
boundaries it must have no dependence on any property that 
differs between the two triangles. Thus slope continuity can 
be achieved only when all slope directions at the boundary 
are free of all properties of the opposite vertex - vertex 1 
in this example. This cannot be achieved by any simple 
weighting of vertices by local coordinates since the 
weighting at the boundary depends on the location of the 
opposing vertex. This is briefly discussed in the Errata to 
McLain (1976). 


A more helpful approach is to look at the work done on 


finite element analysis in engineering. Zienkiewicz (1971) 


“4 


SOA y Whe 7 Coen IA Ae a G 
7 


oy gio t2e%o 0 we tae ye orele | 
Jo hei SAS 1.0 -etae kor ee: aia 34 
cf (was by. ate: wige oat tein bas «Ec 

ue in 2085 rane jon 2% eee r ; a ae “ 


» 
Gabo. ‘anne oat )saerian ait ‘Santa ed a 


a , 
(ey stipe’ (a tih ert BLODCGr4E ano pha ane oe” 
Ace ra Win pak Uns a ao 7 sii hieos isoot ‘gid , ae0>8 é 


OV Vt: we Lj if } fr oO | aot seupay’ chunrtede ‘sunaman 
eile sovaalonh xed ay 242708 ee Rem: enaiaute, 4 ‘3 


fetta 14 (iba Pp a Dieeog sleSfasonsdy tas 


Sie Ve ay orts Tine bY we le find Sim ad ¥ 


Ae pe npy “ 
j cou 


© S Sfbae sive Oba tae ‘tax tot Soap, ar: aaa i i ¢ 


ue ’ 
ogd2 nN?) NOPtones s/ 2T aH TLO eke iaay 


Oe 


Largs e@20ts64 bad hhh dt en Oa: ast 


ted yeaeao, ra VRS LOnopss ‘Hupqes od ‘arth pane 
7 | Se hd i ne me 
aor Va LYE Spo tee, aug? -apipnasas ‘ogi ‘Sia 


ii 
ia BU BAR IFB Sr RIAA EY ka’ ris ‘be 
A ur 


ti aetsaw = veddev @ftadaue aes to ahete 
ie 


at 


gases YWist xd pyvakrng +3 towapy one 


t 4 ark ‘ - bart 
oa: wwe aptent bab. revi py er eis 
of) 0 a mM 
ere ro wine SeSGL aut 0), a ea 
niga 


> Jean 42 b been 
c & ae ge oh ail ~~ cube 


¢ _ 
7 - - : 7 4 


and Irons (1969) are particularly useful, the first of these 
authors giving a relatively straightforward "shape function" 
which is used here. This function can be expressed as 

[6-5] Zp= Disk Lip eji By 
although not used in this form by Zienkiewicz. Here Zp is 


the elevation of the surface at an arbitrary location p and 


the are the coefficients of the planes associated with 


ex 

y) 
the vertices of the triangle, as defined previously. If 
if#j#k #4 i then 


{6.6 ] cjFrj +0. 51; 7 Eke and wij = a ~ Wik 


oN hee 
Where i,j and k refer to the three local coordinates. Ck is 
used here instead of Uae pee in order to preserve the 
symmetry of the operations performed using homogeneous 
coordinates. 

This function possesses the correct slope at data 
points but does not necessarily provide continuity of slope 
across triangle boundaries. Zienkiewicz, however, suggests 
an empirical remedy. If an assumption is made about the 
slope normal to a side at the mid-point of that side, then 
the function may be corrected to conform to this. Thus Sj 
may be defined as the difference between the function slope 
(perpendicular to the boundary) and the average of the 
slopes of the two planes associated with that side, and vj; 
as a function having zero value or slope at sides j and k 
and unit slope at the mid point of side i. Here 

(6.7 ] Vi= 4(V+rj rp Fj SUC A AE | Set (Tj +I,)- 


| 
Addition of Lay 5; to the previous function will provide 
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the necessary slope continuity. Other assumptions can be 
made about the mid-side slopes, producing a variety of 


adequate shape functions. 


6.8 Contour Mapping with the Interpolant 

The surface thus produced is continuous, and is 
continuous in the first derivative, across all triangle 
boundaries. In order to plot the surface, each main triangle 
is divided into N2 equal sub-triangles and Zp calculated at 
each node. Linear interpolation of contour segments within 
each sub-triangle is straightforward since no saddle-point 
problem exists once the basic triangulation has been 
defined. The resolution of the resulting map is a function 
of N and the local data point density. Using this approach a 
contour map was produced (Fig. 20) of the strength of the 
vertical component of the earth's magnetic field, sampled at 
26 sites in Alberta (Bannister, 1977). Fig. 21a Shows a plot 
of sin(x)sin(y) where slopes are known at the data points. 
Fig. 21b is the resulting plot when slopes must be 
estimated, as is usually the case. Apart from boundary 
effects the results agree quite well. These figures 
illustrate the satisfactory performance of the smooth 
interpolating function under real as well as mathematically 
defined conditions. 

It is of interest to compare the results of the 
triangulation method with contour maps produced using the 


grid approach. Fig. 22 shows the bedrock surface of the 


che 


ihe ro Via 7 
; TOT i 
; a GA ee 

/ i ; ; ¢ re Gash 

‘) \ ea y iy it nN 7 7) ia rr, | : 


at ate each scavikien “feat Oley 


o "VrehrTiv « Te — 


| igphoa sett og a . 

sk Hae y a nieiaders NOD) BE reson oak 

fp at Po chie we ore k's ond Sy tage \soaaa 2 
siveseur Abed does ioos Tass ene. sofa, cient 
$e Bossiuefen q° AaB colgunkaendue regia 
RL aetorpse tendo. 10 nob? eL9 ioe | ead 
oy Cau ane FORE babyae teeta be dia Vania 

war ei Shey Pole bak “fh ret 
me csoae, ¥ si tert cu iaeete ait *hO. co 

BS dubue td > puted? sy seated sa box ateh beer 
js 260 WPiette sds 29 0a) \oeeay. Devs teal ‘ome 


oe 


th bal oases yh loca Siyodh on Bee a Cre) 73) 


i. fi i 
tGaq 5 Gaigne ets Ap aa ACCEE vie sat qanm oe . 
2 ho Logs fit Cid 28, TAOS gh 5+ shh aaa. - ) 


iat Setar» meqoike nad 40,L6 er ey 


\ LSI eaerd za os anda 8 As (iiouav aby de 


iat = 


ine 


asinyis eeeit siibs s3tite aoape 2s 
\ 
Asteria (wii Be) a5meug O8nSy concen i 


sigectscusasat a6 Lhew\ap) [saz Gehne ee 


aa 


~ 


“au 
out 2 Geluees one sana 9 papier 
on? pubav Bete hong ro 


eit %» soataue hore eel Mahon sista 


Fig. 


20 (After Bannister, 1977) An example of triangulated 
contour mapping - a magnetic field map. 
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study area produced by the SURFACE II contouring package, 
(Sampson, 1975) from the 78 drill hole data points used in 
this study. Seven of the drill holes used (434E, 725E, 7655, 
809E, 952E, 965E and 774E) were not drilled to bedrock and 
hence the bedrock elevations at those locations were 
estimated visuaily from neighbouring data points. A 50 by 50 
grid was used in association with 1/d2 weighting of a fixed 
humber of nearest neighbours to each grid node. Fig. 23 was 
produced in a similar fashion, but weighting the elevation 
estimate derived from the best fit plane at each data point. 
Fig. 24 shows the same data contoured using the method 
developed here. If these maps are compared with the hand 
drawn map of Fig. 25 which is based on human interpretation 
as well as additional information it will be seen that the 
triangulation method delineates the three major channeis 
much better than the grid approach. Although the hand drawn 
map takes advantage of judgemental data not readily 
quantifiable for computer use, the major features are partly 
resolved by the triangulation of the 78 data points alone, 
There are, however, still several “bridges" where the data 
point distribution imposed a triangulation forcing a ridge 
to cross what is in reality a valley. This is the 
traditional saddle point problem in a different form, and 
these occurrences define regions of inadequate sampling for 
the features to be resolved. However, if a hand generated 
triangulation is defined, based on the additional human 


knowledge (Fig. 26), then the surface may be forced to 
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conform to the predetermined interpretation (Fig. 27). 

This approach is closely related to the idea of surface 
specific lines expressed in Peucker (1972). Surface specific 
points were defined as points that furnish more information 
than merely their coordinates, such as peaks, pits and 
passes. Surface specific lines consist of ridge and valley 
lines, and connect surface specific points, ridges 
connecting peaks and vaileys connecting pits. Passes are 
formed at the intersection of ridge and valley lines. This 
writer knows of no method for analytic definition of ridge 
and valley lines, and hence these must currently be 
identified by the operator and either inserted into the 
program or added to the map at a later stage. 

At least one commercial organization (SCA, 1970) is 
known to use triangulation techniques for contour mapping 
but details of its methods are proprietary and hence remain 
unknown to this writer. No account is taken of the problems 


of surface specific lines. 


-9 Extension into Higher Dimensions 
McLain (1976) suggests the possibility of interpolation 
in higher dimensions but unfcrtunately the author has found 
no interpolating function that guarantees slope continuity 
between domains. Simple planar interpolation however is both 
useful and direct: it is necessary merely to generalize the 


mesh construction for higher dimensions, and a method is 


proposed below. Harbaugh and Merriam (1968) generated 
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hypersurfaces by performing trend surface analysis on crude 
oil gravity values from samples with locations given in 
three dimensions. Similar applications concerned with 
contouring the copper concentration in an orebody, for 
example, could be of geological interest. 

In n-dimensional space a simplex (i.e. the simplest 
polyhedron of that dimension) has nt1 vertices and hence n+t1 
local coordinates, each involving n+t1 order determinants. 
Insertion of new data points into a universal Simplex is 
therefore direct. Optimization of the data structure, 
however, has to be rephrased. An adjacent pair of n- 
dimensioned simplices will have n vertices around their 
common face, and two independent vertices. The alternative 
structure i$ to connect these two independent vertices, 
producing n simplices around this central axis like the 
segments of an eeange This alternative is only possible if 
the simplex pair is convex. 

If this condition is met the optimizatior choice is 
between the simplex pair and the n-simplex "orange", and 
either set may be switched to the other on the basis of the 
maximized minimum height optimization criterion. It is 
interesting to note that coverstg the minimum height is 
one of the few possible criteria inhigher*than-2= 
dimensioned space, since minimizing the diagonal is 
meaningless and maximizing the minimum volume is of little 


benefit when the number of simplices is not the same for 


both alternatives. 
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Appendix 4 contains diagrams of each drill hole in the 
map area. In the centre of each diagram is a reduced copy of 
the electric log for that hole, with broad descriptive 
categories of the material, as well as the location of 
Oxidized zones, in the left part of the central column. The 
calcium carbonate equivalent and bulk chemistry of the 
Ssilttclay fraction are expressed as histograms to the right. 
The far left shows the grain lithology results, and the 
near-left the grain size figures. These last are displayed 
by plotting (from right to left) a dark line for the 10 to 
20 mesh fraction, a gap for the 20 to 40, a line for the 40 
to 60 bate Piowmeanouler gap for the 60 to 140, another line 
for the 140 to 230 mesh size fraction, a gap for the silt 
and a line for the clay fractions. In this fashion much of 
the information needed to distinguish the various till units 
is assembled in one place. The most important additional 
property, lateral continuity, was incorporated by extensive 
correlation between drili-holes, using the properties 
displayed in Appendix 4 as aids in matching till units. 

Examination of the plots in Appendix 4 indicates the 
presence of at least three tills. An upper SPELT Ce yor 
moderate ( 3-6%) carbonate content covering the whole map 
area is underlain in the eastern part of the map area by a 


till of higher ( 5-10%) carbonate content (tL LIS2) so cALthough 
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readily identifiable locally as moderate or high carbonate 
tills, these terms vary in meaning across the map area, all 
tills being higher in carbonate to the east. AS was found in 
Saskatchewan (Christiansen, 1972b), carbonate content was 
the most useful index for till differentiation. Other 
properties, for example Ca, Mg and acid igneous sand grains, 
were closely associated with carbonate content. Additional 
variables were occasionally helpful in indicating till 
boundaries. 

An oxidized zone directly below the upper till was 
observed in one-third of the drill-holes. It thus seems 
possible that sufficient time elapsed between the two 
glaciations for extensive weathering and/or soil development 
to occur. Oxidation of tills by groundwater action does 
occur but in the study area it appears to be sporadic, 
bearing little relation to compositionally-defined till 
boundaries. Such oxidation is usually associated with highly 
permeable zones within the surficial deposits. Additional 
evidence for the presence of a soil includes an oxidized 
zone between two tills, frequently in the absence of highly 
permeable material; the persistence of the zone at one 
stratigraphic position and the presence of a Similar zone, 
interpreted as a weathered surface, in nearby areas of 
Saskatchewan (Christiansen, 1972b). 

Underlying these two tills, but lacking any evidence of 
extensive weathering at its surface, the lower till has a 


low ( 1.5-5%) carbonate content and extends throughout the 
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map area. In the eastern part the electric logs suggest that 
it may be divided into two units (tills 3 and 4), although 
this distinction is not clear on the basis of its chemistry. 
Glacial or preglacial deposits may separate this till from 
the underlying bedrock shale or sandstone. The "Saskatchewan 
gravels and sands" (or “Saskatchewan gravels") are well 
described by Stalker (1968). They include many preglacial 
sand and gravel deposits in Alberta, but the important 
feature is that they are fluvial, lie unconformably on 
bedrock and can be shown to be nonglacial by their lack of 
igneous material derived from the Canadian Shield, their 
abundance of chert derived from the Cordillera and, except 
when located close to the Cordillera, by the absence of 
readily weatherable mineral species. Their presence in 
bedrock valleys is usually considered to confirm their 
preglacial origin. In this study, since interpretation is 
based exclusively on drill hole information, the presence of 
lithologicaily similar material at the base of a bedrock 
valley is considered to suggest a preglacial age for the 
valley . The presence of any Shield material in a clean 
sample of otherwise typically preglacial material would 
strongly suggest the later reworking of the preglacial 
deposits. While it is true that slumpage in the drill hole 
may introduce such material from above, a Similarity of 
grain sizes for both the preglacial and shield material 
would suggest that they were, in fact, mixed prior to their 


being disturbed by drilling. In the plots of Appendix 4 the 
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suggested till units are marked in the right half of the 
central column. The upper till is labelled 1, the middle 
till 2, the upper member of the lower till is labelled 3 and 


the lower member 4, 


J-2 Relation of Analyses to Till Units 

The four till units suggested were distinguished on the 
basis of a wide range of properties, including their 
composition. At least some of the measured variables should 
therefore vary significantly between some of the tills. 
While this may be confirmed by visual inspection on a hole 
by hole basis for some variables, any general statistical 
test, runs into at least two problems with respect to sample 
location. Firstly, sample composition within a till varies 
appreciably and somewhat systematically over the map area. 
Carbonate content, for example, is likely to be lower in the 
western part of the map than in the east, whichever till is 
being considered, and hence it is unreasonable to compare 
directly samples from differing map locations. The second 
problem is that two of the units, the middle till and the 
lower part of the lower till, have only been identified in 
the eastern part of the area, thus biasing their 
compositional estimates if compared with units possessing 
better sample coverage. 

In an attempt to minimize this problem, analysis of 
covariance was performed on each analyzed variable in turn, 


using the ANOVA program available in the SPSS package of 
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statistical programs (Nie et al., 1975). In this a one-way 
analysis of variance is performed to determine if the 
analyzed variable varies significantly between any of the 
four units, but at the same time the sample's location is 
incorporated in the analysis as two covariates (distances 
east and north in metres UTM coordinates). These covariates 
are first applied to explain as much of the variation as 
possible and the analysis of variance performed on the 
residual variance. Linear regression procedures are used on 
the raw data to remove variation due to the covariates, and 
thus more complex variations in composition with location 
are not considered. Nevertheless at least some of the bias 
due to location has been explained. These analyses are shown 
in Appendix 5. The variables were found to vary 
Significantly at the 5% level between the four till units in 
all cases except the 20-40, 40-60, 60-140 mesh and total 
sand fractions. 

If some of the variables differ significantly between 
till units, it is of interest to ask which variables 
differentiate which till units. Table 8 shows the grand mean 
for each variable and the mean for each tili unit, corrected 
for location and expressed as deviations from the grand 
mean. These values are rank crdered with the unit expressed 
as a superscript from 1 to 4, as in Appendix 4. Duncan's 
"new multiple range test" (Steele and Torrie, 1960, Section 
7.5) has been applied to them, and those means found not to 


be significantly different at the 5% level have been 
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Variable Grand Mean Unit Means 

Sand 10-20 3.13 -.183 -.03! 4278 sue 

Sand 20-40 5235 e102) “203! 5017 2419 

Sand 40-60 8474 “1.272 -.113 07% 80° 

Sand 60-140 14.73 -1.49% -.642 -.,063 46! 

Sand 140-230 5.74 -.592 -.35* 018 283 

Total Sand 37.7 -1.372 ~-.74* <-.063 .52? 

Silt 51.37 =1.422 70735 12-66%  _2.39* 
Clay 10.93 ~1.022 -.89*% -.643 290! 

caco3 6.78 ~1.603 -1.27* .138 4,782 
Local 7.07 -3.762 -.163 203* 1.08! 
Acid 59.18 -2.82% -2.223 .98! 4.202 
Basic 4.27 -1.33*% -.203 -.034 .892 

Carbonate 4.04 we62%) (92102 9) 288%) 2 12572 
Quartz 25.44 -2.902 -1.41! 2.27% 2.483 
- 31 SUD Ie i OSC 

ig 1.05 . 4183 -.18* = 0034 502 

Al 2863 etter soy 079 248 

Si 33.11 -2962 = 03% O54 a 27? 

K 2.04 -.033 ~-.03¢ _.01! - 062 

Ca 2.91 | Shen” Ses 2out 2932 

Ti 0.44 -.022. -.01% .023 202° 

Fe 3.82 -.05! 0.02 063 ___ 08° 


of analyses between till units. Column 
for each variable and columns 4 - 
J show the values for each unit, expressed as Se Neca 
from the grand mean. Superscripts indicate the ae 

unit. Underlined means do not differ from each other 


significantly. 


Table 8 Variation 
2 shows the grand mean 
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underlined. Table 9 summarizes which variables may 
satisfactorily distinguish between particular tills in this 
study. Thus a variable plotted in a particular cell has a 
Significantly higher value in the till forming the column 
than in the till forming the row. On the right of Table 9 
are shown those variables whose mean value for a till 
differs significantly from ali the other tills, either by 
being high or low. Those variables in parentheses are second 
highest or lowest in Table 8. It can be seen that ali tiil 
units may be distinguished on the basis of composition, but 
the lower member of the lower till only differs from the 
upper member by having a somewhat coarser sand texture and a 


higher Al content. 
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Clearly the behaviour of some of the variables in Table 
9 is similar to the behaviour of other variables. It would 
be of value for future work to determine which analyzed 
variables group together and hence which ones are redundant. 
Cluster and factor analysis are suitable techniques for 
examining this problen. 

Cluster analysis (Hartigan, 1975, Sneath and Sokal, 
1973) is a method of partitioning an observed Sample 
population into different classes and producing hierarchical 
Classifications of these. In this case the 22 analyzed 
variables themselves were considered the samples to be 


classified. Several different forms of cluster analysis 
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Unit with Higher Value Vaciables Diagnostic for each Unit 
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Table 9 Analyzed variables that differentiate between 
Yort UNnits. AS ar “exanple ; “row? 1", " coltmnr’ 74 “in tthe 
left hand table indicates that Fe is significantly 
pLaiereyin wre ethan Lt unit t. rue crigneenana wtapre 
indicates those variables that differentiate a till 
unit from all others. Variables in parentheses are 
second highest or lowest (see Table 8). 
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available in the CLUSTAN package (Wishart, 1975) were 
performed using the correlation matrix given in Appendix 6 
as the information describing the relations between 
variables. Figs. 28 and 29 show dendrograms produced by 
hierarchical clustering, in which initially each variable 
forms a separate cluster and pinetere are grouped in order 
of greatest Similarity. Similarity between groups may be 
defined in several different ways. Fig. 28 shows the result 
when group Similarity is defined as the lowest similarity 
between any pairs of points, one from each cluster. Fig. 29 
Shows the resuit when the average, rather than the lowest, 
pair similarity is used. Fig. 30 shows the results of using 
the "dendrite" procedure in which the minimum Spanning tree, 
or shortest dendrite connecting all points, (Wishart, 1975, 


p. 67) is broken into the desired number of clusters by 


eliminating the required number of connections. The linkages 


removed are chosen on the basis of minimizing an error sum 
term. A final dendrogram (Fig. 31) is produced by the “mode" 
procedure in which a sphere of a chosen radius is considered 
to enclose each data point, and clusters are defined as all 
those points whose spheres form a common bounding surface. 
The method of introduction of points and updating the sphere 
radius is such that at the start and finish only one cluster 
is present, and it has been suggested that the number of 
clusters present prior to the first fusion of groups 


represents the lowest "natural" level of Classification, in 


this case three clusters. 
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Pig. 28 Dendrogram showing the relationships between 
the 22 analyzed variables. Ordinate values represent 
the product - moment correlation coefficient. 
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HIERARCHICAL CLUSTERING (UNWEIGHTED PAIR GROUP) 


Fig. 29 Dendrogram showing the relationships between the 
22 analyzed variables. Ordinate values represent the 
product - moment correlation coefficient. 
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Fig. 30 Dendrogram showing the relationships between 
the 22 analyzed variables. Ordinate values represent 
the within - groups error sum ot squares. 
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Pig. 31 Dendrogram showing the relationships between the 
22 analyzed variables. val 
product - moment correlation coefficient. 
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The cluster methods illustrated all agree in general on 
these three major groups. Group 1 consists of all sand 
sizes, Carbonate and quartz grains and Si content may 
associate weakly with the sand sizes, but not consistently 
between methods. Group 2 consists of components that could 
be considered "erratic" in origin, possible sub- groups 2a 
containing Na and K and 2b containing acid and basic grains, 
calcium carbonate equivalent, Mg and Ca. The association of 
groups za and 2b with each other is agreed on by all except 
the furthest neighbour method. Group 3 consists of "local" 
constituents -3a consisting of siit percentage and local 
grains, 3b consisting of clay, Al, Ti and Fe. 

Factor analysis (Cattell, 1965a,b, Harman, 1967, 
Kaiser, 1970) is a technique for examining a correlation 
matrix to determine whether its dimensions may be reduced 
without serious less of information. This is achieved by 
creating factors cr components which are linear combinations 
of the observed variables. The first step is to define a set 
of uncorrelated components that explain as much of the 
Variation of the data as possible. The second step is to 
determine if any of these components are of insignificant 
value in explaining the data variation, and eliminating 
them. The third step is to rotate these components until 
they are associated in some fashion with groups of the 
original variables. Appendix 6 contains the results of 
factor analysis of the 22 analyzed variables, using the 


principal compont approach and varimax rotation. 
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Calculations were performed using the SPSS package (Nie et 
al., 1975). Although seven eigenvalues are greater than 
unity, a scree test (Cattell, 1966) on these suggests that a 
maximum of four are meaningful. Component 1 associates 
strongly with the sand size fractions, and negatively with 
the silt fraction. Component 2 associates strongly with 
calcium carbonate equivalent, Ca and Mg, and moderately with 
10-20 mesh sand, acid and basic grains, Na and K. There is 
strong negative association with Si and, to a lesser extent, 
quartz grains. It can be related to erratic components of 
the tili, especially the Devonian carbonates and, perhaps to 
a lesser extent, the igneous Precambrian shield to the 
northeast. Component 3 associates strongly with clay, Al, Ti 
and Fe, suggesting a relationship with the local shales. 
Component 4 associates strongly and positively with Na, K 
and quartz grains but moderately and negatively with Si, 
Suggesting a possible relation to quartz grains and 
feidspars from the Precambrian shield to the northeast. 
Component 5 associates fairly strongly with carbonate and 
quartz grains, Component 6 with the fine sand fractions and 
Component 7 associates positively with acid and basic grains 
but strongly negatively with local grains. 

A comparison with the cluster analysis results suggests 
that components 1,2 and 3 correspond well with cluster 
groups 1,2 and 3, but that there may well be a distinction 


between cluster groups 2a and 2b based on the source of the 


erratic material. 
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Analyzed variables may thus be classified into 
textural, erratic and local groups, implying that in further 
work fewer variables need be examined. In addition, the use 
of factor analysis techniques permits the derivation of 
composite component scores. Thus the Component 2 score for a 
sample may better represent the erratic component of the 
till than an individual variable. While such scores were not 
used in the visual estimation of the till unit contacts 
(Appendix 4) they have been used in the examination of 


compositional variation within a till unit. 


Differentiation 
Examination of the plots of each variable against depth 
Suggests a question: "Which variables can one use in order 
to distinguish different till units within a drill hole?" On 
the basis of cluster and factor analysis the most useful 
variables are those expressing the local/ erratic 
proportions in the tills. Ca, Mg and calcium carbonate 
equivalent are obvious choices that estimate the content of 
erratic material . Clay content, even with the poor accuracy 
of analysis in this study due to long sample storage, 
estimates the local material content fairly weli, but not as 
well as Al or Fee Total sand content, although not 
diagnostic of any particular till units, often is useful to 
mark breaks in depositional continuity. Acid igneous grains 
are also occasionally helpful. However, for identifying till 


units as continuous between drili holes several miles apart, 
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variation in the erratic comronent is the most useful. In 
practical terms, percent total sand and percent calciun 
carbonate equivalent are the variables most commonly used to 
distinguish till layers, and calcium carbonate equivalent 
alone is frequently adequate to correlate these units 
between drill holes. A good example of clearly distinguished 
tills can be seen in drill hole 714E in Appendix 4. It is 
also evident that the bulk chemistry of the till matrix 
(less than 63 microns in this study) exhibits a far greater 
degree of homogeneity within the till units than would be 


expected from examination of the coarser fractions. 


On the basis of the suggested till units and groupings 
of analyzed variables, how should the compositional 
variation in the tills be described? Apart from vertical 
variation in composition between till units there is 
vertical variation within till units at a given location. in 
addition there is the lateral variation of a particular 
till. Boulton and Dent (1974) describe grain size variation 
due to shearing and weathering in the top metre or so of an 
exposed lodgement till, but in the current study the 
variation under consideration is on a larger scale. Appendix 
7 illustrates a property of the lateral variations of the 
analyzed components of the four till units -- it may be very 
large. These maps were produced using the triangulation 


techniques described earlier. If lateral variation is indeed 
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larger than vertical variaticn on a regional basis, then a 
natural question concerns the local area over which lateral 
consistency is observed, since some lateral consistency must 
exist in order to permit the observation of vertical till 
homogeneity. A final question concerns the similarity of 
behaviour, or lack of it, between the different analyzed 


components of the till. 


Hole number 714E was chosen to illustrate the vertical 
Variation in till composition, since three well-sampled tiil 
units were identified in it. The middle till (unit 2 in tae 
analysis of variance tables of Appendix 8) was sampled 19 
times, the upper member of the lower till (unit 3) was 
Sampled 13 times, and the lower member of the lower till 
(unit 4) eight times. Appendix 8 contains the results of a 
one-way analysis of variance performed on each of the 22 
analyzed variables, with the three till units as the 
classifying categories. The main question is whether the 
between-units variation is significantly greater than the 
within-units variation, in which case that particular 
analyzed variable may be used to distinguish between at 
least two of the till units found in this drill hole. 
Appendix 8 shows that this is the case for the 10 to 20 mesh 
sand, all grain counts except the basic igneous grains, 
calcium carbonate equivalent, and all the bulk chemical 


analyses. Calcium carbonate equivalent and calcium and 
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magnesium content showed extremely high levels of 
significance. 

In addition to the analysis of variance tables, 
Appendix 8 shows the mean, range and 95% confidence interval 
of the mean for each of the till units. From these it can be 
determined which tills may be distinguished using a 
particular analysis. The middle till (unit 2) may be 
distinguished quite well from both the lower tills by the 10 
to 20 mesh sand, local, acid igneous and quartz grains, 
Calcium carbonate equivalent and all buik chemistry 
analyses. The upper and lower members of the lower till 
(units 3 and 4) may be distinguished by calcium carbonate 
equivalent, Mg, Si, Ca and Ti concentrations. In addition, 
erratic components tend to increase in concentration, and 
local components decrease, with decreasing age of the till 
unit. This is a variation between units, rather than a 
simple variation with depth, since for many analyses the 
between-units variance wasS much greater than that within 
units. Hole 714E shows compositional differences between 
tills better than most of the drill holes. It does however, 
illustrate the high possible homogeneity of the composition 


of a till unit throughout its thickness, especially in the 


fine fraction. 
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In contrast, the lateral homogeneity of tiil 
composition is far less clear, and the composition maps of 
Appendix 7 may show appreciable fluctuation in composition 
over the study area. This suggests that till forming 
processes were local in effect but stable throughout the 
period of deposition of at least the basal part of a till 
unit, although distinction of basal from ablation till has 
not been attempted in this project. Sufficient information 
is not available to suggest a depositional mechanism that, 
while only producing the same effects over short distances 
laterally, maintained these effects during most of the 
depositional period of the till. It appears unlikely that 
such a process would remain locally stable over extended 
periods of time while the flow conditions of the ice sheet 
changed in response to climatic and other conditions. It 
appears much more likely that the actual deposition of till 
only occurred as a single relatively short-lived phenomenon 
in the life-cycle of a continental ice-sheet. Under these 
conditions the buik of the till unit + up to 100 ft (30m) in 
the case of hole 714E, and greater in some others - would 
have been deposited under one, perhaps very short-lived, set 
of local flow conditions. While deposition may have occurred 
more than once, under different local conditions, there is 
little evidence for this in the upper and middle till units, 
although this is a possible explanation for the weak but 


sometimes measurable differences between the upper and lower 
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members of the lower till. 


If depositional processes were involved that remained 
constant over only limited lateral extent it is important to 
examine the extent of those possibly local processes. An 
attempt was made to investigate this by using a suitable 
analysis of variance design. In addition, since some till 
constituents were described as being erratic components, it 
was necessary to determine whether they or any other 
analyzed variables could be considered to vary significantly 


between different parts of the study area. 


7.9 Sampling and Statistical Model 

In order to compare local and regional lateral 
compositional variation it was necessary to select several 
distinct areas of the map each containing several drill 
holes. In addition, if between-units variation was to be 
examined these drill hoies should penetrate several tills, 
and if vertical variation is cf interest, several samples 
should be available from each unit. Unfortunately, since 
computational methods become extremely complex for suitable 
analysis of variance designs if some categories contain no 
samples, the choice of test data was very limited. 
Consequently three areas were selected, a western one using 


holes 435E, 436E, 437E and 758E. A central area containing 


holes 456E, 498E and 709E; and an eastern area using holes 
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457E, 458E and 765E. Since the middle till and the lower 
member of the lower till appear to be absent from the west 
of the area, oniy the upper till and the upper member of the 
lower till were used. This was perhaps unfortunate since the 
middie till was compositionally the most distinctive. 
Finally, due to the low number of samples in some drill 
holes only two samples were used from each unit in each 
drili hole, selected from near the top and bottom of the 
unit. 

The analysis of variance design used was a model 
described by Winer (1962, p. 334, Table 7.3-10). The model 
would have been a simple factorial design were it not for 
the fact that a particular drill hole clearly cannot be 
repeated in each of the three areas. Hence drill holes are 
"nested within" areas and are treated as replicates. 
However, the unit and position within the unit of each 
sample must still be accounted for and are "crossed" with 
each other and with the drill holes, since a sample may be 
found for each combination of drill hole, unit and sample 
position. This model is most easily defined as a single 
factor model (area) with replicates (drill holes) and 
repeated measures. These repeated measures therefore consist 
of the four combinations of unit and sample position. It is 
a property of this model cand of the data) that interaction 
terms involving drill holes cannot be tested and must be 
pooled with other terms to form the relevant error sums of 


Squares. Thus the sum of squares due to area isitested 


TTS 


BAT 


sewgi) od 4 Bas ates cinta ae 

n 

So oe BaP wots * tecanee &2. & ad oF Hatta - 
gue to abudeer wena add: sranahs gee 


Ske ongee 2 Akt toRe  DIsdse4g eb wee 


reoubeeby deo Pr, atowors 
ae ae arose ot elie Va sedate, obec 
cans lL Ts - soph as -). Deseo on 
edt. HONMES TO. ia? got ete ean ante 
| | ae a 
fohoa & waw +eop ch heed sun a | 
igiom pT .(RPsbel eldet SOekiae Nae Sie 
“eet wn. eh Stee © i alt | Luknouans idan 
{ -wsabo videaloveled Wbbte : rshuokra ag. & ae 4 

az4 2alod Li Feb, wo0S8- sees seat ‘pis os , bie 
,eaeeokige? 28 auth wy fe dig: aepean 
inee to, dite as ChAT se ILI TAGg” bas: tka 
thaw: PReeae eal ie ap Bele, 2e3 aay 0 sttne d 
at sine ii eal ‘sg ponte >. 2816n (ike . in) ied, fad 
sig mee. Oehe’ Sia. Shon “EL agh Zo ae anne ee 2 
eignka 6 jaa featinb Yibiave san ak, depos bce 
hak denhor rae: aby perenblgee Hai ere 


$2ieneo oysie rsd? det was fadeaqad saad? val 


| Bie 41! a ORTLAO ston Dae TEM: 10) ti 
f>igosreseL ted? waesvais to Bie) pba ee so) 
tH sau@ine dos63% od gouges Sibu: Lipo gaketoad 

ta save Joins MmaAVeLer @HtimhOR os ened 2agison ds tu 


bovae? ai.sots ot oul goteppe 20) o 


against a pooled error sum of squares consisting of the sun 
of squares due to drill holes plus the interaction between 
drill holes and areas. This interaction has no conceptual 
meaning since drill holes are nested within areas. Similarly 
the significance cf units, and the interaction between areas 
and units, are tested against the interaction between units 
and drill holes plus the interaction between units, drill 
holes and areas. Sample positicn and (area x position) 
interaction are tested against the (drill hole x position) 
plus (drili hole x position x area) interactions. Finally, 
the (unit x position) and (area x unit x position) sums of 
Squares may be tested against the pooled (unit x position x 
drili hole) plus (unit x position x drill hole x area) sums 
of squares. 

Appendix 9 shows the analysis of variance results for 
each of the 22 variables and Table 10 summarizes the 
significance of the relevant F, or variance ratio, tests. A 
particular concern is the relation between local and 
regional variation for each component of the till. Those 
analyses differing significantly between regions of the map 
as opposed to within regions were acid igneous grains, 
calcium carbonate equivalent and Ca content at the 5% level 
and Na at the 10% level, with Mg approaching 10% 
significance. Thus, although the variance analysis was 
performed on fewer samples than would have been desirable, 
it appears that only erratic constituents of the till vary 


more on a regional basis than a local one within a specified 
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till unit. Appendix 7 shows the lateral variation in each 
till unit for grain size, grain count, calcium carbonate 
equivalent and bulk chemistry. 

The till components that do vary significantly between 
regions of the map are those that correspond strongly with 
Component 2 in the factor analysis described in Section 7.3. 
It was suggested that this component represented erratic 
material introduced into the map area by the moving ice. As 
can be seen from the maps of Appendix 7 these same 
components, despite any local variations, decrease generally 
towards the southwest. Besides maps of calcium carbonate 
equivalent for each till unit, Figs. 32 to 35 contain maps 
for each till unit of the Component 2 scores for all 
samples, calculated from the factor-score coefficient matrix 
of Appendix 6. The score for each sample is obtained by 
normalizing each analysis value (by subtracting the mean and 
dividing by the standard deviation for that analysis), 
multiplying it by the Component 2 score coefficient for that 
analysis, and then summing these results over all analyses. 
In addition Figs. 36 to 39 show maps of the ratio of the 
concentrations of erratic cations (Na, Mg, K, Ca) over local 
Sations (Al, Si, Ti, and Fe) that_are probably derived 
mainly from the local sandstones and shales. Caliciun 
carbonate equivalent, Component 2 scores and erratic ratios 
tend to behave in a Similar fashion for a particular tiil 
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variation in carbonate content as a function of grain size 

and distance of transportation from the source. They found 

that with increased distance the original relatively coarse 
material was subject to more crushing, so that the maximum 

concentration of erratic material was found in successively 
finer size fractions until a stable "terminal grade" was 


reached. Thus conclusions based on one size-fraction do not 


readily relate to those based on a different fraction. Since 


both calcium carbonate equivalent and bulk chemistry are 
measured on the silt plus clay fraction in this study the 
regional variation is cleariy due to dilution by local 

material in the down-ice direction rather than by further 


grinding action. This holds true for all till sheets. 


Due to poor sample distribution available for the 
analysis of variance of Appendix 9 it was not possible to 
Compare the middle till, with its distinctive high erratic 
components against the other units. The two more extensive 
tills that were compared were fairly similar in their 
composition, since only the 40 to 60 mesh sand and the clay 
content differed significantly at the 5% level between ali 
units over the map area. The sand figure is difficult to 
interpret, but the clay figure reinforces the visual 
observation that clay values often changed appreciably 
between units. Significant variation between the top and 


bottom samples of a till was only found for Na at the 5% 
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level, and for 20 to 40 mesh sand, silt and Mg at the 10% 
level. Considering the small number of samples and the low 
number of significantly varying analyses, it seems unlikely 
that a clear general change in composition of a till can be 
expected with depth. This agrees with visual impressions 
from the drili hole plots in Appendix 4 and the analysis of 
Appendix 8. 

Various interaction effects were also found 
Significant. The significant interaction of unit and sample 
position for 40 to 60 mesh sand, and carbonate grains, 
suggests that some units may be stratified with respect to 
these properties, but the lack of corroborative evidence 
makes this inconclusive, although the interaction between 
area, unit and depth found significant for sand and clay 
measurements suggests some textural stratification in some 
units in some areas. The interaction of area and unit is of 
interest, since it was found significant for all size 
analyses except 140 to 230 mesh sand, and clay. This 
strongly suggests that some areas of some units were 
markedly different in overall grain size from others, and 
that this effect, although affecting each till unit 
independently, was extensive enough to apply toeaiset wot 
drill holes rather than just one. Examination of individual 
values shows that the upper till tends to be particularly 
coarse in the westernmost of the three areas examined, and 
the upper member of the lower till coarser in the 


easternmost area. Thus the texture of individual tills may 
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be fairly homogeneous not merely throughout their depth, but 
laterally for several miles, This interaction suggests that 
local depositional conditions responsible for depositing 
those till units that are compositionally fairly homogeneous 
may have an extent roughly comparable to the areas chosen 
for the analysis of variance design (10-20 km). Examination 
of the grain size maps in Appendix 7, as well as the 
analysis of variance results, suggests that these mediun- 
scale variations in composition may be associated with the 
topographic features encountered by the moving glacial ice. 
Gillberg (1965, 1967a,b) did extensive work on the 
relationship of the erratic component of till to the 
underlying topography. He clearly showed an exponential 
decrease of erratic material with increasing distance from 
the source material, primarily due to dilution with local 
material. He also demonstrated the appreciable variation in 
local ice direction in the lower part of an ice sheet (and 
hence in local till composition) in response to sub-ice 
topography, even though the upper, major portion of the ice 
sheet moved in the regional direction. He concluded that if 
the topographic slope was such as to impede the ice flow the 
erratic material wculd be deposited rapidly, whereas with a 
down-ice slope, erratic material would be transported 
greater distances. This corresponded to up-ice dips 
producing zones of deposition of both the local and the 
erratic material being transported, and down-ice dips 


producing zones of erosion of local material. In the lower 
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member of the lower till in the study area for this work the 
erratic component tends to be lower in the Helina and 
Beverly Channels and higher on the St. Lina slope, 
suggesting that, according to Giliberg, the former would be 
zones of easy flow and the latter a zone of resistance. The 
Kikino channel did not exist at that time. In the upper 
member of the lower tili the erratics tend to be low in 
parts of the Beverly and Kikino Channels, suggesting fairly 
free ice flow in these, but much of the Helina channel shows 
littie influence on the till composition, as much of it had 
previously been filled in. 

The middle till varies widely in its erratic content 
and displays little pattern, and the upper till, although 
decreasing fairly systematically in erratic content towards 
the southwest, tends to have low values on the high ground 
and higher values in the residual low areas. No obvious 
explanation is available for this, although the erratic 
behaviour of the middle till may be due to its termination 
in the map area and consequently varying modes of deposition 
could well have been in effect at various times and places. 

Thus, although medium-scale lateral variation in till 
composition may well reflect local depositional conditions 
controlled by the effect of the underlying topography on the 
continental ice sheet, more detailed sampling would be 


necessary to confirm this. 
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Two Upper Cretaceous bedrock formations directly 
underlie the surficial deposits in the map area (Carrigy, 
1972). These are the Lea Park and Belly River Formations 
(Fig. 40). They dip gently to the southwest, with the 
younger Belly River only present in the southwest quarter of 
the map area. The Lea Park Formation consists of massive, 
marine, silty, dark grey shale and argillaceous siltstone 
with some fine-grained sandstone lenses. The Belly River 
Formation consists of non-marine, fine- to medium- grained, 
grey argillaceous bentonitic sandstone, grey clayey 
siltstone and grey or carbonaceous shale (Yoon and Vander 


Pluym, 1974). 


Westgate (1968) cited McConnell (1886) as recognizing 
the presence of elevated, commonly gravel-capped plateaux in 
the Plains of Western Canada The implication was that the 
plateau surfaces once formed the level at which rivers 
flowed, and uplift and subsequent degradation lowered the 
surrounding country some 600 m. Because drainage was to the 
east, material clearly of Cordilleran origin was deposited 
throughout the plains. These formed the Saskatchewan Sands 
and Gravels which are well described by Stalker (1968). 


Although their lithology indicates this origin, rock and 
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Fig. 40 {From Yoon and Vander Pluym, 1974) Bedrock 
geology of the Edmonton - Lac La Biche - Cold Lake 
area. 
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mineral species resistant to weathering and erosion 
predominate. The immediate lccal sources of these gravels in 
the Edmonton area were the late Tertiary conglomerates that 
covered former extensive preglacial uplands (Rutherford, 
1937). Much of this material was deposited in the preglacial 
valleys as the uplands were denuded. Stalker (1968) provided 
criteria for their recognition, including the requirement 
that they lie between the bedrock and the first glacial 
deposits. These gravels in the buried valleys form a major 


groundwater source in the map area. 


Three major buried channels have been located in the 
test area by previous workers - the Beverly, Helina and 
Kikino Channels. (Yoon and Vander Pluym, 1974). The Helina 
Channel is extremely broad, with low slopes to the valley 
walls, its shape suggesting a mature major river valley. The 
Beverly Channel is less broad and has steeper bedrock walls 
(see Fig. 25). Preglacial material has been identified in 
the laboratory from the Helina Channel (drill hole 833E) and 
the Beverly Channel (706E and 831E). Field identification of 
preglacial gravels without laboratory confirmation occurred 
in holes 808E and 809E in the Helina Channel. The Kikino 
Channel is fairly narrow, steep-walled and cuts across the 
heck between the other two channels. No preglacial material 
has been found in the Kikino Channel, although a field 


identification in hole 730E was rejected in the laboratory 
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as being contaminated with Shield material. 

To the east the merged Helina and Beverly Channels pass 
through Coid Lake and join a known bedrock channel that 
drains east (Christiansen and Whittaker, 1974). The known 
Vermillion Channel which enters the southeast corner of the 
map area also joins this river system in Saskatchewan. The 
smaller Vegreville Channel joins the Beverly just east of 
the Kikino Channel. The Beverly Channel is followed by the 
current North Saskatchewan upstream from the map area. 

The Beverly Channel presumably contained the preglacial 
North Saskatchewan River (Westgate 1969, p.133). The Helina 
Channel probably contained the preglacial Athabasca River 
and possibly even the headwaters of the Peace River which 
May at one stage have drained east through Lesser Slave 
Lake, Calling Lake and Lac La Biche. (see Fig. 1 and 
McCrossan and Glaister, 1966) The Kikino Channel is thought 
to be of an early ice-marginal origin, draining water from 
the ice front. ann the North Saskatchewan northwards through 
Lac La Biche. The region around Lac La Biche possesses some 
very thick (greater than 100m) surficial deposits of gravel 
that is either of glacial or reworked preglacial origin. 
This valley deepening and gravel accumulation may be due to 
the Kikino ice marginal channel intersecting the preglacial 
Athabasca or Athabasca / Peace river system. Further 
drilling would have to be done to the north and west to 
resoive this problem. Surficial topography suggests that the 


preglacial Athabasca may have entered the map area from the 
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WNW, via Fawcett Lake and Calling Lake. The Kikino ice 
marginal channel may have continued NNW via Wandering River 
and Pelican Portage, associated with the breach in the 
mountain range now occupied by the Athabasca River flowing 
north towards Fort McMurray (see Fig. 1). 

Recent work by Yoon et al. (1977) indicates the 
presence of an additional channel, the Primrose Channel, 
running just to the north of the study area between Imperial 
Mills (just to the north of drill hole 958E in Fig. 2) and 
the northern edge of Cold Lake (near hole 765E). Preglacial 
Sands and graveis have apparently been identified in it. The 
preglacial Athabasca River may therefore have flowed east 
from near Lac La Biche to Cold Lake by either the Primrose 
Or ne Helina Channel at any given time. Further work is 
needed before confident descriptions may be made of the 


preglacial drainage in the area. 
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The sequence of maps of Figs. 25, 27 and 41 to 59 
Summarizes the glacial history of the area. Fig. 25 shows 
the current bedrock surface drawn by hand, utilizing not 
only information derived from the 78 drill-holes in Appendix 
4 but data from other drill-holes and bedrock maps from 
other sources (Carlson 1967, Yoon and Vander Pluym, 1974). 
Fig. 27 is a computer drawn map of the same surface produced 
from the 78 drill-holes alone, as are the rest of the maps. 
The Helina and Beverly channels drained east through Cold 
Lake prior to glaciation, but the Kikino Channel in the west 
had not yet been cut. Fig. 41 shows the current thickness of 
materials deposited underneath the first glacial till whose 
bottom contact is shown in Fig. 42 and thickness in Fig. 43. 
This material includes preglacial sands and gravels in the 
Helina Channel as well as material deposited during the 
glacial advance. It can be seen from Fig. 43 that the Helina 
Channel was partially filled with till at this time, and the 
Beverly channel may have been treated similarly. The ice 
depositing this lower member of the lower till seems to have 
reached most of the way up the St. Lina slope (at the 
junction of the Beverly and Helina Channels) and stopped. 
Since this would have blocked the drainage of the preglacial 


North Saskatchewan it was forced to cut the narrow ice- 
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marginal Kikino Channel in order to drain northwest through 
Lac La Biche. The preglacial Athabasca river would have been 
treated similarly and the water of both dammed rivers, 
together with the ice run-off would have been channeled 
northwards towards Fort McMurray as described earlier. 
Although it is usually considered that ice-marginal 
channels are supplied predominantly by ice run-off, the 
Kikino channel does not appear to be present to the south of 
the Beverly channel. Although ice-shoved bedrock is common 
in the area, the bedrock topography map (Yoon and Vander 
Pluym 1974, Fig. 6) does not indicate even a blocked route 
for a southern extension of the Kikino channel, strongly 
suggesting that this channel is in fact due to redirected 
river drainage. This idea was also made by Rutherford (1941) 
to explain the generai drainage pattern of the Saskatchewan 
River system, where most channels show an abrupt bend to the 
southeast suggesting that drainage, being blocked by ice, 
was forced to flow along or roughly parallel to the ice 
front. Ultimately the channels turn again to the northeast, 
indicating that the regional slope was still in that 
direction. In the same paper the author showed the 
horthwestern extension of the "Altamont" moraine associated 
with the Coteau in Saskatchewan passing along the more 
recently discovered Kikino channel. Another mcorainal system, 
the "Viking", is shown passing in a north-south direction a 
few miles further to the west. This suggests that one, if 


not two, regional glacial advances stopped near the western 
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portion of the map area. 


2 Deposition of the Upper Member of the Lower Till 
Following the deposition of the lower member of the 
lower till the ice probably retreated to the northeast. On 
the basis of the estimated topography of the top contact of 

the lower member of the lower till (Fig. 44), the Helina 
Channel seems to have remained blocked. The Athabasca 
continued to drain north, but it is uncertain whether the 
North Saskatchewan resumed its old course as suggested by 
Fig. 44, or continued to follow the Kikino Channel as 
Suggested by the bottom contact of the upper member (Fig. 
46). Fig. 45 shows the thickness of non-till material below 
the upper member and above the lower member in the eastern 
part, where the lower member ice did not advance. It is 
frequently impossible to tell, in the absence of a 
particular tiil in any one drill hole, whether the non-till 
material was deposited before, after or contemporaneously 
with the missing till. Consequently Fig. 45 and other maps 
of this type assume the material to be as young in age as 
the tills present will permit. In the case of Fig. 45 the 
deposits between the upper and lower members of the lower 
till therefore include in the west some preglacial gravels 
as well as some presumed glacio-fluvial material deposited 
in front of the ice. 

Fig. 47 shows the thickness of upper member of the 


lower till deposited. It can be seen to cover the whole map 


o - ri : : } “i : - 7 ) 


Os tecellraou o> otf bianeel eae sot ois 
26 ‘goedeos, gox ony to qin tpoded bad 
| end foie i CPW ry 2p) ds iia - 
yqaiited tn out cbesosiOe boadsiter noua 
F wr 
‘ ade yeueeuy Rha PP Bey ee ili atl oid 
ed Regeeppre to aatogn Boag i. Sone ae 
ee Toitvand: Yaa LM ‘ad? oral Leek. Oey ue O 
ee 2) 
ape wi vodwem aqow ait do +ontaos soto oy ate 
vohad (be breeae 2 L300 i vasrnatolt he oie mene) 
ifes2ee Sida ol. T9deon nsinbil out gvodE Bai 
eth iaomerbs SPOR Bet le By schist nena 
Rhian bal gi Py ae stint 
(fhe gon ed? Sedsvadw iis lod hay shi ane. ate pee 
; Ae a 
= y bia ah: Ai nies ies , an nay ie@ xoeee ite 
ae ' ea ; 
es xe bas CH ieRl “/ se Ue ROD 
io Beal ae BiIuYy 66 aq) OF Preerg y) 
aa7 eb. O28 is PREO =? ‘= ~.veeiog 
raeul aa Le a ee RON OM daa: a3qau 4 
. Si HZE neus ba albus ch tod ween 
é o » : VATS 
I 7 ee! b 
ay OF paren ae Lebte ten nbd sobaed p bamuined 
cn os ij es - 
i i hoy aud 4% ‘aedawe OH; att” ‘tia | 
ee ea “ei ne wai anid pve co ine c 5 
mip 
haraiy. Uy . . Aik "ie 
‘ ie 


L i - a i? Ae a mon ry) 


ee 
u? 
i 


A tn Np 


_——— gs 


146 


3 
3tun TTT} HburdqTaepun sToaexyetpoeuut 
Ss} tsodep PpezepTt Tos uooun 
JO 3808F UT SSeuyosTyL Gh ‘*hHhta 


bs it a es, 
' | | 


mG ee De a we 
oy a yn y | 
ga pa rire veut 


7 


147 


ey 


3 


*€ 9TuUn TTT32 JO eseqg 
JO YO UT VOTPeASOTA Oh °Hhtg 


RE 


= 


ou 7 15a 


i 3 


Bre" ike BIS4SS 


Ve 
i 
iS 
ia 
é 


arma 


gerne 3 7 


Tae 


or 
bench. ea ooh a 


149 


area with the exception of the Whitefish Lake bedrock high. 
Thicknesses of greater than 300 feet (100m) are found 
filling in the recently-formed Kikino Channel. This great 
till thickening suggests that the higher ground on the 
western boundary of the map area again formed an impediment 
to further ice movement. That this is the case is supported 
by the fact that some or much of the upper member till on 
the west side of the Kikino channel is composed of disturbed 
bedrock. This is the case for drill-holes 9525, 953E, 956E 


and 957E. 


9.3 Bedrock Thrusts 

St. Onge (1972) noted the hummocky moraine associated 
with the high ground just to the west of the region now 
associated with the Kikino channel. He remarked ({ p. 93) 
"Ridges on both sides of Lake Amisk and south of Kikino are 
also associated with it. The ridges south of Kikino 
represent a large push moraine constructed by southwesterly 
flowing ice." While the current work suggest that the ridges 
were formed prior to the deposition of the uppermost till, 
both sets of ridges are associated with disturbed bedrock - 
the Lake Amisk ridges are near drill-holes 952E and 953E 
while the Kikino ridges are near drill-holes 956EF and 957E. 
These morainal features may well be associated with the 
Suggested Viking moraine of Rutherford (1941). 


Of particular note are the sampled drill-holes 434F and 


828E together with the other drili-holes in the cross- 
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section of Fig. 48. On the basis of drill holes on either 
Side of the cross-section (Yoon and Vander Pluym 1974, Fig. 
6) the Beverly channel should intersect the cross-section 
near drill hole 434E with a valley floor below 1800ft 
(550m). However, the drill holes (which are numbered in 
chronological order) prior to number 828E had identified a 
bedrock plateau at 1850ft. between holes 435F and 436E, and 
a Surprising bedrock high of approximately 2100ft. elevation 
where the channel was expected. The subsequent drilling of 
hole 828E to depths well below the first occurrence of 
bedrock material initially indicated that permeable material 
occurred below 130ft. of shale. This suggested that the 
missing valiey fioor had been found at this location, with a 
rather implausible bedrock high immediately to the north. 
Subsequent examination suggested that the bedrock high was 
in fact an ice-shoved bedrock slab of shale overlain by 
sandstone. The Lea Park Shale - Belly River Sandstone 
contact subcrops in this area (Fig. 40), and holes 759E and 
796E to the south terminate in the sandstone while holes 
735E, 754E and others to the north terminate in the shale. 
For this reason the first sandstone / shale sequence in 
holes 828E and 434E is interpreted as an upward-thrust ice- 
shoved bedrock slab. The second sandstone is interpreted as 
Belly River Sandstone in place, underlain by Lea Park Shale. 
Hole 795E also terminates in the Belly River Sandstone. It 
is thought that this bedrock slab was thrust across the 


Beverly Channel by the ice advance that deposited the upper 
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member of the lower till. This interpretation is supported 
by the presence of bedrock thrust ridges further north in 
the Kikino Channel (St. Onge, 1972) and also by extensive 
disturbed bedrock found in drill holes 952E, 953E, 956E and 
957E. 

The high ground formed by the bedrock thrust extends 
nearly to Wahstac and appears to have formed a permanent dam 
preventing the North Saskatchewan from following the Beverly 
Channel any further (Figs. 25 and 49). The river, thus 
dammed, was forced to cut an outlet to the south, eroding 
over 100 feet (30m) of bedrock in the process. It is not 
clear to what extent the. upper member of the lower till 
extends to the west of the Vilna-Wahstao barrier. 

While the separation of the lower till into upper and 
lower members may be uncertain on the basis of composition 
alone it is guite clear that some action, prior to the 
dominant upper member deposition and subsequent to the 
initial disruption of the preglacial drainage, must have 
caused the Kikino Channel to be formed, and a lower member 
ice advance as estimated from electric logs reached 
approximately the right location. An alternative, if 
similar, interpretation is that the upper member glaciation 
paused for an extended period of time just to the east of 
the current Kikino Channel position and the ice marginal 
channel was formed, only to be buried immediately the ice 


readvanced. 
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Fig. 50 shows the elevation of the top of the lower 
till. There is no evidence of major drainage in the area, 
although local streams may have flowed east to Cold Lake, 
Neither extensive leaching nor clear zones of oxidized iron 
associated with this till are found in the Sand River area, 
Fig. 51 shows the inter-till deposits directly below the 
middle till, consisting mainly of valley-fill deposits near 
Cold Lake, and Fig. 52 shows the pre-middle till surface. As 
seen from Fig. 53, the thickness of the middle till is 
greatest ee old Beverly Channel and at the top of the 
St. Lina slope, to the west of which it was not found in the 
drill-holes. This thickening, along with its topographic 
expression as shown in Fig. 54, suggests the formation of an 
end-moraine, possibly the Altamont of Rutherford (1941). 
Fig. 55 shows the inter-till deposits below the upper till - 
the four patches of thick deposits are clearly deposits 
whose true age is older but can not be so identified due to 
the absence of the required till layer at that location. Fig 
96 shows the topography at the base of the upper till, and 
most of the weathered zones located are associated with this 
Surface, indicating an extended ice-free period. There is no 


obvious drainage pattern. 
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Fig. 57 shows the thickness of the upper till, which 
extends over the whole area, and also continues further to 
the west. The predominant effect is of further deposition of 
till to reduce irregularities in the underlying terrain, 
together with a tendency for the till to thicken at the 
western edge of the map area where the ice met the higher 


ground. 


& Present Topography 

Fig. 58 shows the elevation of the top surface of the 
upper till. It suggests that little major drainage passes 
through the area, although as shcewn on Fig. 59 and Fig. 2, 
(the generalized map of the present topography), a broad 
vestigial valley is still found following the Helina Channel 
between Lac La Biche and Cold Lake. The current rivers, 
primarily the Amisk, Beaver, and Sand Rivers, are relatively 
Small with very localized drainage basins. A minor low 
region, occupied by Whitefish, Goodfish and Garver Lakes 
remains as evidence of the Kikino Channel and its effect on 
ice movement, but little evidence of the Beverly Channel 
remains. Indeed, a conspicuous neck of high ground still 
crosses the channel between Vilna and Wahstao, evidence of 
the magnitude of the Vilna bedrock thrust. Just upstream 
from this barrier, St. Onge (1972) noted a system of 
meltwater channels that may have drained the North 


Saskatchewan River to the north at one stage in the retreat 
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of the most recent glaciation. Ellwood (1961) suggested a 
deglaciation sequence for the Vermillion area to the south 
that takes into account the meltwater channels in that area, 
as well as the late-stage remobilization of stagnant ice 
that produced the south-east trending flutings that are such 
a conspicuous feature on air-photographs. At an even later 
stage of deglaciation, Kocaoglu (1975) noted glacio- 
lacustrine deposits just below the junction of the Beverly 
and Kikino channels, and fluvio-glacial deposits Slightly 


further upstream. 


9.7 Historical Summary 
The glacial history of the area may be summarized with 
reference to the fence diagram of Fig. 49. 
a7 
The preglacial North Saskatchewan and Athabasca 
Rivers flowed east through the Beverly and Helina 
channels, joining in the eastern part of the map 
area. 
2/ 
Glacial ice advanced into the map area from the 
north-east, depositing till unit 4 in the eastern 
part of the map area and disrupting the regional 
drainage. The North Saskatchewan River, being dammed, 
was diverted north, cutting the Kikino Channel. The 


Athabasca River was probably diverted north, out of 


the map area, at this time. 
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3/ 
The ice retreated to the east, and while the North 
Saskatchewan may have resumed its old course, 
abandoning the Kikino Channel, the Athabasca did not 
re-occupy the Helina Channel within the map area. 

7 
After what may have been a relatively short time 
interval the ice readvanced into the area and 
deposited till unit 3. In the process of overriding 
the newly formed Kikino Channel bedrock thrusting 
occurred, completely blocking the Beverly Channel 
upstream from its junction with the Kikino Channel, 
forcing the North Saskatchewan to cut its present 
channel to the south. It is not known how far the ice 
extended to the west of the Kikino Channel. 

5/ 
Following the retreat of the ice an indeterminate 
period of time elapsed before the arrival of the next 
ice sheet. No major drainage passed through the area. 

6/ 
The new ice advance deposited the distinctive high- 
carbonate till unit 2 over the eastern half of the 
area, but apparently did not advance any further to 
the west. 

if 
Following the retreat of this ice sheet an extensive 


period of time elapsed prior to the advance of the 
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final ice sheet, as indicated by the extensive 
weathered zone formed in many places on the land 
surface of that time (Section 7.1). 

8/ 
The final ice advance from the east extended past the 
western boundary of the map area and deposited till 
unit 1. After the ice retreated this till formed much 


of the present topographic surface. 


2.8 Correlation of Surficial Units 

Christiansen (1968a,b, 1972a,b) and Whitaker and 
Christiansen (1972) classified the surficial deposits of 
central Saskatchewan on the basis of drill-hole information. 
The Empress, Sutherland and Saskatoon Groups established by 
them are described below. 

The Empress Group is defined as those deposits between 
the bedrock surface and the lcwest till, including both 
preglacial and glacial material. The Sutherland Group 
consists of the lowest major till and is distinguished from 
the overlying tills by possessing more clay and less sand, a 
generally lower electric log resistivity and a lower 
carbonate content, especially of dolomite. The Sutherland 
Group can sometimes be separated into an upper member with 
moderate and a lower member with very low carbonate content. 
Partial leaching of carbonates is sometimes found in the top 
40 feet of the Sutherland Group. 


The overlying Saskatoon Group is composed of the lower 
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Floral till and the upper Battleford till. These are not 
usually separable on the basis of carbonate content, but the 
Floral is frequently oxidized towards the top, indicating 
fairly extensive weathering, and is more compact and jointed 
(with oxide staining in the joints) than the overlying 
Battleford till which is more friable, soft and unjointed,. 

Table 11 shows the tentative correlation between 
Christiansen's work in Saskatchewan and the results of this 
study. On the basis of a lower carbonate content, lower 
electric-log resistivity, frequently lower sand values and 
the occasional separation into two units, the lower till of 
the Sand River area is correlated with Christiansen's 
Sutherland Group, although leaching of carbonates at the top 
is less obvious in the Sand River area. On the basis of the 
occurrence of an oxidized weathered zone directly below the 
upper till, the middle and upper tills of the Sand River 
area are correlated with the Floral and Battleford tills 
respectively. It should be noted that frequently in 
Christiansen's published cross-sections the Floral till does 
in fact have a rather higher carbonate content than the 
Battleford. In agreement also is the somewhat higher Mg 
content in the upper till of the Sand River area by 
comparison with the lower till. 

With respect to carbonate values it should be noted 
that Christiansen's cross-sections tend to be parallel to 
ice-front positions and therefore approximately equidistant 


from the source of the erratic carbonate. Hence he does not 
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specifically consider the decrease in carbonate content with 
distance from the source. In the Sand River area however, 
the main channels are roughly parallel to the ice flow 
direction and the effect of diluting the erratic carbonate 
with local shales is marked. This may well be due to 
relatively high erosion of local material due to the 


regional uphill slope facing the advancing ice. 
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It has been shown in this work that it is possible to 
decipher the till lithostratigraphy and broad glacial 
history of an area from drill-hole information. The tools 
required are primarily the traditional ones of driller'ts and 
electric logs. The visual examination of cuttings or side- 
wall samples is not enough, and some compositional data are 
required. Different tills probably differ in composition in 
moSt areas because succeeding glaciations mask or expose 
different portions of bedrock, and the earlier tills 


themselves may be incorporated in the younger ones. 


10:2 Choice of Compositional Measurements 

Analysis of various compositional rarameters indicates 
that many tills are remarkably homogeneous throughout their 
thickness, notwithstanding the occasional textural variation 
or sand lens that is apparent from the electric log. This is 
especially true of the silt and clay fraction in the present 
Study. Lateral variation is less consistent, being affected 
by local sub-ice conditions and topography, although 
Systematic dilution of erratic materials is evident away 
from their source. Compositional parameters should therefore 
be chosen from the fine rather than the coarse fraction of 
the till, since more individual grains will usually be 


analyzed in a procedure performed on a fine fraction. Total 
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sand content is frequently useful in separating till units 
within the drill-hole, as is the electric resistivity log. 
Compositional parameters in the fine material heavily 
reflect the ratio of erratic to local material, and in this 
Study, as in others, calcium carbonate equivalent has proved 
to be the most concise summary of this, as well as being one 


of the easiest analyses to perforn. 


On the basis of the plots in Appendix 4, a few 
empirical suggestions may be made on sample collection in 
the field. Several samples should be taken from each 
suspected till unit since a single sample may not be 
representative. When plotted in the fashion of Appendix 4, 
four samples wili probably visually indicate the average 
composition of the till unit. The sampling decision would 
normally be made in the field on the basis of the 
previously-run electric log. Breaks in the general trend of 
the resistivity curve frequently denote a new lithologic 
unit and form a good basis for field sampling decisions. 
Samples every 5 feet are recommended for thin units, and 
every 20 or more feet for very thick tills. Deeper tills 
Should be sampled at least as carefully as more shallow 


ones. 
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4 Regional Stratigraphy 

Four possible tills have been identified in the study 
area. Several of them may be regional in extent, as has been 
Shown in Saskatchewan. It is Suggested that Christiansen's 
(1968a,b, 1972a,b) nomenclature for surficial deposits 
(Empress, Sutherland and Saskatoon Groups; Floral and 
Battleford Formations) be used whenever the interpretation 


seems justified. 


Till stratigraphy has some benefit for estimating the 
hydrologic potential at various depths in an area. Sand or 
gravel layers are less likely to occur within a till than 
between tills, and if present are less likely to be 
laterally continuous. The zone between the lowest till and 
bedrock has a higher potential for water supply than the 
Overlying tills, especially if associated with a bedrock 
valley. Inter-till sands are also fairly extensive, but 
usually tend to be less clean. Aquifers are most likely to 
be found in depressions in the surface on which they were 
formed - whether in a bedrock valley or a superimposed till 
Surface. A basic understanding of the glacial stratigraphy 
and history should therefore facilitate the estimation of 


groundwater potential. 
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Techniques for bulk carbonate analysis, bulk chemistry 
of the fine fraction of till samples, and the computer 
display of contour maps have been developed for this study. 
in each case these developments appear to be tools the use 
of which extends beyond this work. The equipment developed 
is a step in the direction of automated laboratory carbonate 
analysis in both surficial geology and soil science. Bulk 
chemical analysis on fused till samples is currently 
dependent on the availability of an electron microprobe, but 
Should soon be a viable technique on the more common, and 
less expensive, electron microscope when fitted with a 
Suitable detector for energy-dispersive analysis. In terms 
of the rate of analysis of multiple elements and samples the 
method compares favourably with other available techniques. 

Finally, the developments in automated contour mapping 
described here provide both facilities and economies not 
previously available. The approach outlined has several 
advantages in terms of fiexibility, economy and the 
properties of the surface representation produced. Future 
enhancements could easily include slope calculations, 
automatic hill shading, storage of multiple surfaces (for 
example, geological contacts) and volumetric and areal 
calculations. In contrast with many methods the ideas used 
(1.e. breaking the region into triangles, obtaining the 
Slopes at each data point and bending triangular plates to 


conform to these) are easy to comprehend, if not always easy 
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to implement, and hence the casual user should not receive 
any unpleasant surprises in the form of outrageous costs or 


surface shapes along with his desired Map. 
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research in the Cold Lake oil sands area. Canada - 
Venezuela oil sands symposium 77 proceedings. (Awaiting 
publication). Branch Rept., Alberta Environment, 
Edmonton. 


Zienkiewicz, 0.C., 1971. The Finite Element Method in 
Engineering Science. McGraw Hill, London. 
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Crush air-dried sample. 

(Due to lengthy storage the sidewall samples are 
usually extremely hard. Consequently coarse crushing should 
be done with a jaw crusher followed by hand crushing with a 
rolling pin to avoid breaking rock fragments. Continue hand 
crushing until little new material passes through a #40 
Sieve.) 

Sieve on #40 sieve (using "Sonic Sifter" for one minute 
at setting 6). Weigh and record the +40(A) and -40 (B) 
fractions. 

Wash +40 fraction on a #40 sieve and oven-dry. 

Sieve on #10, #20 and #40 sieves using Sonic Sifter. 
Weigh and record the +10(C), +20(D) and +40(E) fractions. 


Store -40 and total +40 fractions. 


Sieve +40 fraction on +10 and +20 sieves. Use -10, +20 
fractions. 

If necessary, partition sample by quartering until at 
least 200 grains are left in the pile. Categorize 200 grains 
as local shale or sandstone, acid igneous, basic igneous, 
Carbonate or undifferentiated quartz. Record number of 


grains in each category, recombine and store +40 fraction. 
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Sieve -40 fraction using #230 sieve until 5 to 10 gm of 
-230 fraction are obtained. Weigh and record. Place this 
sample in approximately 250 ml of Calgon solution (40 gm 
Calgon per litre distilled water), agitate and leave 
overnight to soak. 12 to 24 samples may be prepared at one 
time. 

Distilled water should be kept at room temperature, 

Wash sample into soil dispersion mixer with a smail 
amount of distilled water. Mix at low speed for two minutes 
and high speed for one. 

Wash sample back into jar with distilled water, make up 
to 400ml, agitate thoroughly and place on solid surface. 
Record starting time. 

Calculate settling time for 2 micron fraction to fall 4 
cm at room temperature using Stokes Law. Assuming a 
particle density of 2.65, T in minutes = 10¢k2 where 
K=0.01365 at 20°C. 

At end of settling time insert the pipette to the 4 cm 
mark below the surface of the liquid. Siphons off 25 ml, 
place in a weighed bottle and oven-dry. 

Record dry weight of the aliquot after subtracting the 
weight of the bottle and the weight of Calgon in a blank 
aliquot. 

Wash remaining sample in a #230 sieve. Oven dry +230 
Geaction. 


Sieve on #60, #140 and #230, sieves. Weigh and record 
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the -40 +60(H), -40 #¢140(I) and -40 #230(J) fractions. 

Since it is necessary to use small samples and to 
extract some of the -230 fraction for further analysis, data 
reduction is as follows: 

- a correction factor due to loss of -230 fraction, 

k = 1+ (A-E)/B is calculated. 

- the fine fractions (H,I,J) are multiplied Dy K. 

- the silt + clay weight = k (A-J). 

- the clay weight = volume of sample suspension x clay 


aliquot weight / aliquot volume. 


Carbonate Analysis 

Power on the carbonate apparatus, and allow 10 minutes 
for the detector to stabilize. The recorder should be set at 
approximately 5 volts full scale and 8 minutes per inch. 

Obtain 1 gm of sample by quartering the preserved -230 
fraction and place it in the rinsed acrylic flask. Put 10 ml 
of 50% HCl in the reservoir, add a magnetic stirring bar, 
screw down the lid and allow the internal pressure to 
stabilize for a minute. 

Invert flask to mix acid and sample. Place flask 
upright on magnetic stirrer set to about half speed. 

Remove flask and unscrew lid when chart curve indicates 
no further pressure rise. Measure and record pressure rise 
due to the sample above the pressure in the uninverted 
sealed flask. 


Run 5 calibrations of 200 mg oven-dryed (105°9c) calciun 
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carbonate reagent at least once a day, or if marked 
barometric changes occur. 

Calculate % calcium carbonate equivalent for each 
Sample assuming the average pressure change of a batch of 


calibrations to be due to 20% of calcium carbonate. 


Obtain approximately 2 gm of preserved -230 fraction by 
quartering. Place in clean vial with 2 parts by volume of 
lithium tetraborate flux. Mix and pour into a clean 
crucible. Compact sample by tapping and invert onto graphite 
holder. 

Place sample in image furnace, turn furnace up to 60 
volts and allow sample to melt for 10 minutes. Remove 
Sample, quench and break off tip of sample. Place in epoxy 
mount made by drilling 12 3/16" holes through a 1" disc of 
Araldite with a polyethylene ("Saran Wrap") film on one 
Side. Record sample position. Fill mount with epoxy and 
leave 2 days to cure. 

Grind on paper lap using #240 and #400 sandpaper. 

Re-impregnate surface with 3M epoxy. Place in a vacuun 
for 10 minutes. Heat in oven at 110° C for 3 hours or 
Overnight at 72°-75°A4cC, 

Grind briefly on #600 sandpaper, then wash it carefully 
with soap and water. 

Impregnate mount with Geonite and wash with warm water. 


Grind on paper using 6,3 and 1 micron diamond paste for 
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1, 1 and 3 hours respectively. Wash and re-impregnate with 
Geonite between paste changes. 

Clean mount thoroughly, polish with .05 micron alumina 
for one minute, wash with warm water and soap. 


Carbon coat mounts. 


installation) 

Place mounts in the ARL-AMX electron microprobe 
equipped with an Ortec lithium-drifted silica energy 
dispersive detector. Standards to be used are Kakanui 
kaersutite for Na, Mg, Al and Ca, quartz for Si, Hohenfels 
Sanidine for K and ilmenite for Ti and Fe. These must be 
analyzed at the beginning of each run, and a calibration 
willemite standard should be processed at the beginning and 
end of the day's run. Instrument settings are 15 kV 
operating voltage, 300 nanoamps beam current and 8 microns 
scanning diameter for the beam. Beam current and probe 
current readings should be taken every hour with the Faraday 
cage. 

Eight regions of each glass should be counted for 50 
seconds each, for a total of 400 seconds per sample. Data is 
collected on the Texas Instruments 733 cassette tape 
terminal, and the saved data transmitted to the Amdahl 
470/V6 computer for processing through the EDATA energy 


dispersive data-reduction program. 
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APPENDIX 2 


Calcium carbonate equivalence and grain size analyses 
for all analyzed samples. 


The absolute values of the clay figures are considered 


to be unreliabie due to disaggregation problems, 
although usefuvi for comparison with neighbouring 
samples. . 


Samples are identified by drill hole number and depth 
in feet. 
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APPENDIX 3 


Grain lithology and bulk chemical analyses for all 


analyzed samples. 


Samples are identified by drill hole number and depth 


in feet. : 
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APPENDIX 5 


Analysis of covariance for each analyzed variable, 


comparing between - units and within - units variance. 
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APPENDIX 6 


Principal components factor analysis of the an 
analyzed variables, showing the correlation matrix, 
the pre- and  post- rotation factor matrices, the 
transformation matrix, final communalities for each 
variable, eigenvalues for each factor and the 


factscore coefficients. 
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APPENDIX 7 


Contour maps showing the lateral variation of sand, 


Clay, grain lithology, calcium carbonate equivalent 
and bulk chemistry in each of the four till units. A 


heavy black line marks any till unit boundaries. 
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APPENDIX 8 


Between-units versus within-units analysis of variance 
for each analyzed variable, with statistics describing 
the variation within each unit. 
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APPENDIX 9 


Analysis of variance for each analyzed variable, 
_ comparing variation between areas, Pils (Anuss and 
vertical positions within tali-anits, 


The western test area contains drill holes 435E, 436E, 
$37e@ and “/58E, the central area contains holes 456F, 
498E and 709E, while the eastern area contains holes. 
4578, 458E and J65E. 
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